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TITULO 

COMPOSICI6N COHERENTE DE SEflALES POR CORRECCiGN FOCAL 
PROGRESIVA. 



5 SECTOR DE LA TtCNICA 
Tecnologfa fisica 

Sistema electronico para Ecografia clfnica y Ensayos No Destructives. 
Otras aplicaciones: radar, sonar y seismologia 



1 0 ESTADO DE LA itCNICA 

En diversos ambitos de aplicacion, un conjunto de N transductores o array, 
recibe senales de una o multiples fuentes ubicadas a cierta distancia y 
orientaci6n relativa al array. Esta situaci6n es comun en radioastronomi'a, 
radar, sonar, sismologia y, particularmente, en el ambito de la tecnologfa 
15 ultrasonica basada en arrays (ecografia clfnica o industrial, en este caso 
aplicada a los Ensayos No Destructives), donde N se encuentra, tfpicamente, 
en el rango 32 a 256 o m£s. 

En estas situaciones, cada uno de los N elementos del array recibe la senal 
generada por cada fuente con cierto retraso temporal respecto a los restantes 

20 elementos. La funcion de un conformador de haz es obtener una unica senal a 
partir de las N recibidas en la que los retrasos temporales para cada posible 
distancia en una orientacion determinada se compensen, aumentando la 
capacidad de discriminacion del sistema (mayor resolucion, relacion senal/ruido 
y contraste). Repitiendo este proceso en una diversidad de direcciones 

25 mediante barridos iineales o angulares, se obtiene informacion suficiente para 
poder componer una imagen 2D o 3D de las fuentes de senal. En el caso 
particular de los ultrasonidos, las N sefiales recibidas son ecos producidos por 
discontinuidades en el medio que ha sido insonificado por un emisor, que 
puede ser el mismo array receptor (o un subconjunto del mismo) u otro 

30 transductor (array o monoelemento). 
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En recepcion, es deseable compensar los retardos en todos los puntos o focos 
a lo largo de la li'nea de propagation en una direction detemninada. Esto 
requiere modificar continuamente los retardos en el tiempo, para adaptarlos a 
la llegada de los ecos desde distancias dinamicamente variables a cada 
5 elemento. Esta tecnica de focalizacion dinamica, proporciona la mayor 
resolucion, contraste y relation senal/ruido. 

Siendo tecnicamente imposible modificar en forma continua los retardos, es 
esencial que su discretization se realice con una excelente resoluci6n temporal 
para cada elemento del array y para cada distancia, tipicamente del orden de 
10 1/16 a 1/64 el periodo fundamental de la serial. 

Este hecho ha motivado la busqueda de metodos para la focalizacion dinamica 
en recepcion sencillos y eficientes que, a lo largo de los anos, se han ido 
adaptando a las disponibilidades tecnologicas. Asi, se han descrito tecnicas 
basadas en dispositivos CCD accionados por relojes cuyo periodo determina el 
15 retardo [Macowski, "Ultrasonic electronic lens with reduced delay range", US. 
Pat Nr. 4.058.003, Nov. 1977] o en Imeas de retardo conmutables [Maslak, US. 
Pat. Nr. 4.550.607, "Phased array acoustic imaging system", Nov. 1985]. Estos 
m6todos introducen artefactos de conmutaci6n, son voluminosos y 
proporcionan una resolucion limitada estando, practicamente, en desuso. 

20 Otra alternativa, asimismo analogica, se basa en la utilization de mezcladores 
y desplazadores de fase que realizan los retrasos temporales en el dominio de 
la frecuencia [Maslak et al., "Dynamic focused linear phased array acoustic 
imaging system", US. Pat 4.699.009, Oct 1987]. Uno de los problemas 
asociados a este tipo de tecnicas es la dificultad de su adaptation a 

25 transductores de distintas frecuencias. 

Por otra parte, se ha propuesto la utilization de rotadores de fase, basados en 
el algoritmo Cordic [R. Andraka, "A survey of CORDIC algorithms for FPGAs", 
Proc. 1998 ACM/SIGDA VI Infl Symposium on FPGAs, Monterrey, CA, pp. 191- 
200, Feb. 1998]), y su aplicacidn a la focalizacibn dinamica se encuentra en [O 1 
30 Donnell et at, "Method and apparatus for digital phased array imaging", US Pat. 
Nr. 4.983.970, Ene. 1991] y en [Engeler et al. "Dynamic phase focus for 
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coherent imaging beam formation", US Pat. 5.111.695, May 1992]. El hardware 
involucrado es significativo. 

La disponibilidad de tecnologia VLSI permitio nuevas propuestas basadas en 
retardar las senales una vez digitaiizadas. Cada elemento del array se asocia a 
5 un canal de adquisicion y procesamiento (focalizacion) independiente. Las 
salidas de los N canales se suman para obtener la serial resultante (suma 
coherente). Para operar en tiempo real, cada canal dispone de dos espacios de 
memoria: uno para datos, donde se almacenan las muestras adquiridas, y otro 
de focalizacion que establece los retardos a aplicar. 

10 En un principio, la resolucion de los retardos se limito al propio periodo de 
muestreo de las senales, determinado por el criterio de Nyquist, alrededor de 
1/3 a 1/5 del periodo de la serial, utilizando registros de desplazamiento de 
longitud variable [Welles II et al., "Architecture for ultrasonic imaging", US Pat. 
4.796.236, Ene 1989]. La baja resolucion temporal produce Idbulos de 

15 cuantificacion de retardos que reducen el rango dinamico, contraste y relation 
senal/ruido de la imagen [D. K. Peterson, G. S. Kino, "Real-time Digital Image 
Resconstruction: A Description of Imaging Hardware and an Analysis of 
Quantization Errors", IEEE Trans. Sonics Ultrasonics, 31, pp. 337-351, Jul. 
1984]. 

20 La option de utilizar una mayor frecuencia de muestreo no es aconsejable por 
el mayor coste de los componentes, consumo energetico y volumen de datos a 
almacenar y procesar. Para solventar este problema se han propuesto diversas 
tecnicas que, basicamente, pueden agruparse en dos: 

1) Realizaci6n de un retardo grueso con una resoluci6n temporal equivalente al 
25 periodo de muestreo, y otro fino mediante interpolacidn de la serial recibida 
entre dos muestras consecutivas, con una resolucion, tfpicamente, entre 1/4 y 
1/16 del periodo de muestreo. 

A este grupo pertenecen las patentes: [Odell, "Digital beamforming and filtering 
circuit", U.S. Pat. Nr. 5,268,877, Die. 1993], [D.Lipschutz, "Time multiplexed 
30 digital ultrasound beamformerUS.Pat. 5,469,851, Nov. 1995], [Beaudin et al. 
"Delay generator for phased array ultrasound beamformer M , US. Pat. 5.522.391, 
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Jun. 96], [Wright et al. "Method and apparatus for receive beamforming 
system", US Pat. 6.042.547, Mar. 2000] y [M.H. Bae, "Ultrasonic signal focusing 
method and apparatus for ultrasonic imaging system", US Pat. 6.231 .51 1 , May. 
2001]. 

5 Esta tecnica, con diversas variantes para la realization ffsica de los retardos 
grueso y fino, incluida la de rotation de fase mediante algoritmos CORDIC, 
requiere invariablemente una etapa de interpolation y un selector de muestras 
interpoladas por cada canal, con la consiguiente complejidad y errores 
asociados. 

10 2) Adquisicion directa de las senales en los instantes que corresponden a la 
llegada de la serial de cada foco a cada elemento. Para ello, se controla la fase 
del reloj de muestreo de cada canal con una resolution temporal equivalents 
entre 1/16 y 1/64 el periodo fundamental de la senal. Este metodo hace 
innecesaria la etapa de interpolation, eliminando la electronica y las fuentes de 

15 error correspondientes. Las muestras se pueden adquirir sobre una memoria 
FIFO, de la que se extraen para realizar la suma coherente [Chang et 
al.,"Bandwidth sampling technique for digital focusing in array imaging 
systems", US. PAT. 5,581,036, Die. 1996]. 

El mayor problema asociado a esta tecnica es la realization de los generadores 
20 de reloj de muestreo de fase controlada con la suficiente precision. En 
[O'Donnell et al., "Method and apparatus for fully digital beam formation in a 
phased array coherent imaging system", US Pat. Nr. 4.809.184, Feb. 1989] y 
[Corl, "Digital ultrasound system with dynamic focus", U.S. Pat. 4,974,211, Nov. 
1990], se generan m relojes derivados de un reloj maestro de frecuencia m 
25 veces superior a la de muestreo, cada uno de ellos con un desfase 2n/m; un 
selector del reloj por cada canal, accionado por el contenido de la memoria de 
focalizacion, pemnite asociar en cada instante uno de los relojes al de 
muestreo. La distribution de multiples relojes y su selection con coherencia 
temporal es un problema, especialmente para las mayores resoluciones 
30 temporales. 

Otra posibilidad es almacenar el reloj de muestreo de cada canal en una tabla 
como una secuencia de 1s y 0s, como se ha sugerido en [J.H. Kim, T. Song, 
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S.B. Park, "Pipelined Sampled-Delay focusing in ultrasound imaging 
systems, Ultrasonic Imaging 9, pp.75-91, 1987], pero no es practica por la 
cantidad de memoria necesaria para almacenar todos los relojes de muestreo 
de todos los canales para multiples Ifneas de barrido. 

5 Siguiendo otra alternativa, se calculan en tiempo real los instantes de muestreo 
mediante circuitos especializados asociados a cada canal. Esta tecnica se 
describe en [Park et al. t "Real time digital reception focusing method and 
apparatus adopting the same", US Pat. 5.669.384, Sep. 1997], [Bae, " Focusing 
delay calculation method for real-time digital focusing and apparatus adopting 
10 the same" US Pat. 5,386.881, Feb. 1998] y [Petrofsky, "Method and apparatus 
for distributed focus control with slope tracking", US Pat. 5.724.972, Mar. 1998]. 
Todas, con diferentes variantes, se basan en el algoritmo del punto medio 
descrito en [Bae et al., An Efficient Real Time Focusing Delay Calculation in 
Ultrasonic Imaging Systems, Ultrasonic Imaging, 16, pp. 231-248, 1994]. 

15 Aunque se trata de una tecnica interesante por sus posibilidades de integration 
en un ASIC, presents diversos inconvenientes entre los que destacan: la 
necesidad de programar un gran numero de parametros para cada li'nea de 
barrido y cada canal, la complejidad del hardware para obtener una buena 
resolution, menor eficiencia de utilization de area de silicio frente al requerido 

20 por sistemas basados en memoria por la baja regularidad del hardware y, 
finalmente, requiere un compromiso entre obtener una alta velocidad de calculo 
o una reduccidn del rango de aplicacidn. 

El metodo y sistema electronico objeto de la presente invention pertenece a 
este grupo y resuelve los problemas descritos mediante una tecnica mixta, que 
25 combina una memoria de focalizacion de reducido tamano, en la que solo se 
requiere una fraction de bit por cada foco y canal, con el calculo en tiempo real 
del instante de muestreo con elevada resolution mediante circuitos sencillos y 
regulares de facil integracidn VLSI. 



30 
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DESCRIPCI6N DE LA INVENCI6N 
INTRODUCClON: 

Se describe un conformador de haz que realiza la composition coherente de 
las senales procedentes de una direction de propagation determinada, 
5 recibidas por un conjunto de N transductores receptores o array. La 
metodologia y aparatos descritos son aplicables en campos como imagen 
ultrasonica (ecografia y Ensayos No Destructives), radar, sonar, aciistica, 
seismologia y otros. 

El conformador de haz aqui descrito opera con arrays cuyos elementos pueden 
10 estar distribuidos arbitrariamente. En particular, realiza la focalizacidn dinamica 
con deflexion del haz en acimut y elevation utilizando arrays bidimensionales. 

Basa su operacibn en muestrear la serial recibida por cada elemento del array 
en los instantes que corresponden a los de llegada de la serial desde multiples 
focos situados a lo largo de una direction de propagation dada. Dentro de los 

15 limites establecidos por los principios que guian la presente invention, el 
instante de adquisici6n de la senal procedente de cada foco y recibida por cada 
elemento del array se determina con un error absoluto, respecto al valor 
exacto, inferior a medio periodo de un reloj maestro. La frecuencia del reloj 
maestro puede elegirse para reducir la incidencia negativa de los lobulos de 

20 cuantificaci6n de los retardos en la relaci6n senal/ruido, contraste y rango 
dinamico de la serial resultante. 

La frecuencia de muestreo es independiente de la de la serial, siempre que se 
verifique el criterio de Nyquist, tanto en radiofrecuencia como en banda base, y 
de la frecuencia del reloj maestro. Con tecnicas de sobremuestreo (conversion 
25 AL), la frecuencia de digitization puede ser la mas baja compatible con la 
relation senal/ruido deseada. En particular, puede utilizarse el reloj maestro 
modificado segun el metodo descrito para realizar directamente el muestreo 
AE. 

Al reducir la frecuencia de muestreo a la mas baja posible, compatible con el 
30 ancho de banda de las senales y/o relation senal/ruido, se reducen el coste, 
consumo energ6tico y volumen de datos a procesar. En particular, con el 
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metodo que constituye una parte esencial de la presente invention, solo se 
adquieren las muestras necesarias para realizar la suma coherente, sin que 
sea necesario un proceso de interpolacibn o la adquisicibn de informacibn 
redundante. 

5 Ademas, el metodo de conformation que se describe y que constituye otra 
parte esencial de la presente invention, posee unos bajos requisitos de 
memoria de focalizacion, que pueden llegar a ser de una fraction de bit por 
cada muestra adquirida en cada canal. Esta alta eficiencia facilita que el 
metodo pueda aplicarse convenientemente a una diversidad de modos de 
10 operation, como son el barrido lineal o sectorial, con aperturas densas o 
dispersas, y otras aplicaciones en las que se requieren multiples descriptores 
de la focalizacion dinamica para cada configuration de la apertura activa del 
array. 

Dichas caracteristicas posibilitan la integracibn conjunta de la memoria de 
15 focalizacion, la de adquisicion de datos, la logica de generation de los relojes 
de muestreo y la de otras funciones auxiliares como apodizacion, apertura 
dinamica y control, de multiples canales, en un unico dispositivo VLSI o 
estendar configurable tipo FPGA, con la consiguiente reducci6n de coste sin 
perdidas de calidad. 

20 

Descripcion del metodo. 

Por sencillez, la exposition se realiza para un array unidimensional de N 
elementos, no necesariamente equiespaciados, considerando una aplicacion 
de ultrasonidos utilizando la tecnica de pulso-eco en un medio homogeneo con 
25 velocidad de propagation c. El mgtodo que se describe a continuacibn se 
denominard gen6ricamente tecnica de Correccidn Focal Progresiva o CFP, que 
describe la forma en que opera. 

La Figura 1 muestra un diagrama de la geometria general con los 
condicionantes anteriores. El array, de N elementos, numerados 1, 2, 3, /c, 
30 AM y A/, se situa sobre el eje de las x, con el origen de coordenadas 
arbitrariamente definido, por comodidad, en el centra del array. Se desea hacer 
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una adquisicion de serial en la direction angular 9 referente al eje z. Para ello, 
se definen multiples focos, a intervalos AR en la direccitin de propagacibn 
(llegada de la serial). El m6todo propuesto combina las senales procedentes de 
cada foco y recibidas por los diferentes elementos del array de forma que se 
5 obtiene una salida que representa la serial que sen'a proporcionada por un 
transductor del tamano del array focalizado en todos y cada uno de los focos. A 
este proceso se le denomina composicidn coherente de las senales recibidas y, 
al resultado, serial focalizada dinamicamente con deflexion. 

Los intervalos AR se eligen convenientemente en funci6n del periodo de un 
1 0 reloj maestro Tx=1/fx y de la velocidad de propagation c como: 

&R=vcT x /2 (1) 

donde v es un numero entero que expresa el numero de periodos del reloj 
maestro que hay entre dos focos consecutivos en pulso-eco. Por otra parte, los 
intervalos AR deben ser lo suficientemente pequenos como para que se 

15 verifique el criterio de Nyquist: si T s es el intervalo temporal entre muestras, y 
Tr es el periodo fundamental de la senal, el criterio de Nyquist se verifica, 
tipicamente, para Ts < Tp/3, considerando el ancho de banda de las senales 
habituales y la adquisicion en radiofrecuencia. Preferentemente, el intervalo 
temporal entre muestras es un multiplo del periodo del reloj maestro, esto es, 

20 T s = sTx, con s entero. 

En general, puede adquirirse un numero m de muestras entre dos focos 
consecutivos, eligiendose preferentemente el valor v en funcion del intervalo 
entre muestras T" s = sT x para que v = ms t con my s enteros. En una 
realizacion particular en la que m=1, cada muestra coincide con un foco, de 

25 modo que todas las muestras se focalizan dinamicamente. Cuando m>1, 
estrictamente solo se realiza la focalizacion dinamica en los focos, aunque la 
profundidad de foco permite, en general, suponer bien focalizadas las muestras 
intermedias. Una variante del metodo, denominada tecnica de Correction Focal 
Progresiva Variable o CFPV, varia dinamicamente el numero de muestras 

30 adquiridas entre dos focos, reforzando la idea de mantener una buena 
focalizacion dinamica en todas las muestras adquiridas. 
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Por otra parte y, como es bien conocido, el valor \x = Tr/T x o relacion entre el 
periodo de la senal recibida y la resolution temporal del sistema, afecta al 
rango dinamico y contraste de las imageries construidas con multiples haces 
focalizados. En general, se considera aceptable un valor \i > 16 (ver, por 
5 ejemplo, [D. K. Peterson, G. S. Kino, "Real-time Digital Image Resconstruction: 
A Description of Imaging Hardware and an Analysis of Quantization Errors", 
IEEE Trans. Sonics Ultrasonics, 31, pp. 337-351, Jul. 1984]). 

Es importante que la resolucibn temporal determinada por 7> sea 
independiente del periodo de muestreo T s que se elige en base al criterio de 

10 Nyquist y, por tan to, se relaciona con el periodo de las senates recibidas, T R . 
Cuando el muestreo se realiza en radiofrecuencia se elige, en general, T s » 
Tr/3o T s « Tr/4 para reducir el flujo de datos a procesar y el consumo 
energetico, esto es, la frecuencia de muestreo mSs baja compatible con el 
criterio de Nyquist. En tal caso, con s > 6 o s > 4, respectivamente, se tiene 

15 que |i > 16. Una particularidad de las tecnicas CFP y CFPV es que permiten 
operar a la mayor parte de la electronica a la frecuencia de muestreo, 
reserv£ndose la mas elevada del reloj maestro para muy pocos elementos que 
realizan el ajuste del instante de muestreo. 

Para un foco F, situado en la coordenada polar gen6rica (R, Q) y, para el 
20 elemento k del array situado en la coordenada cartesiana (x k , 0), por el teorema 
del coseno (Figura 1 ): 



donde L k (R,Q) representa la longitud del camino recorrido por la senal desde el 
foco F situado en (R,Q) al elemento k. El tiempo de vuelo de un pulso 
25 ultrasonico desde su emision arbitrariamente elegido en el origen de 
coordenadas al foco F y de este al elemento k, es: 



En el caso en que el origen de tiempos se situe en una coordenada diferente 
(por ejemplo, cuando el espacio inspeccionado sea iluminado desde un emisor 




(2) 




(3) 
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diferente del propio transductor receptor), el primer termino de la Ecuacion (3) 
debera ser modificado para representar el tiempo de vuelo desde el emisor al 
foco F. 

En general, se denomina T ki al tiempo invertido por la senal correspondiente al 
foco Fj en alcanzar el elemento k del array. Esta notacion generica se refiere, 
tanto a la geometna indicada en la Figura 1, como para otras situaciones en las 
que los elementos del array se distribuyen de forma arbitraria, el medio de 
propagacion no es homogeneo o los focos se situan en otras posiciones. No 
obstante, para facilitar la exposition, se mantiene la geometna de la Figura 1, 
el medio homogeneo con velocidad de propagacion c y las posiciones de los 
focos se alinean formando un angulo 0 respecto al array situandose, de 
momenta, a intervalos regulares AR a partir de una distancia inicial R 0 en la que 
se coloca el primer foco F 0 . 

La variation en el instante de llegada al elemento k de senates procedentes de 
dos focos consecutivos F h1 y F,- situados a distancias R r AR y Rj del origen de 
coordenadas, respectivamente, esta dada por la diferencia: 

AT h (R. ,Q) = T k (R.,Q)-T k (R. -AR,Q) (4) 
que, con AR suficientemente pequeno, es: 



dTAR,Q) 

Ar 7 .(^. 3 e) = M— £ 

' dR 



(5) 



R = R. 



20 Derivando la Ecuacion (3) respecto aRy sustituyendo el valor de AR dado por 
la Ecuacion (1): 



'J; 



R. -x.senQ 
i k 



vT 



2 + x 2 k - 2R.x k senQ 



X 



(6) 



Esta es una funci6n monbtona creciente entre R=Q y R-* oo e independiente de 
la velocidad de propagacibn, con los siguientes limites: 
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lim I 




1- 



senQ 




>0 



(7) 



V 



lim \ATj=vT 



(8) 



donde, por claridad de notation, AT k = AT k (R,B). El valor xi</\x k \ de la Ecuacion 
(7) representa el signo de x kt mientras que 9 interviene con su signo. En 
5 cualquier caso, el menor valor posible para AT k es cero, independientemente 
de la position del elemento k. 

De este modo, los intervalos entre los instantes de muestreo que corresponden 
a dos focos consecutivos estan limitados por los valores expresados en las 
Ecuaciones (7) y (8). En particular, para el caso 9=0 (no se realiza deflexi6n), la 
10 funcion crece monotonamente entre vTy/2 y vfxcuando la distancia a la que se 
situan los focos varfa de 0 a °o. Cuando se reciben las senales con deflexion 
del haz, el rango de variation es mayor y, en un caso extremo de deflexion a 
90°, entre 0yv7 x . 

En una posible realization, las posiciones de los focos quedan definidas por 
15 valores comprendidos en el intervalo [0, v], que representan los valores 
incrementales dados por la Ecuacion (6) en periodos de reloj maestro. En este 
sentido, si el numero entero v puede expresarse con b bits, cada retardo 
incremental para cada foco y para cada elemento puede representarse con 
esta misma anchura, independientemente de la position ocupada por el 
20 elemento del array, de la distancia a la que se situe el foco y del angulo de 
deflexion. El error absoluto con el que se determina el instante de muestreo es 
inferior a medio periodo de reloj maestro al redondear al entero m£s proximo el 
resultado de ATid(R t Q)/Tx dado por la Ecuacidn (6). El inconveniente de esta 
opci6n es que requiere un elevado numero de bits para expresar los retardos 
25 incrementales para cada foco y elemento, con las consiguientes demandas de 
capacidad de memoria. 
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Inversamente, si se dispone de b bits para representar cada retardo 
incremental, el rango de valores posibles puede elegirse entre dos situaciones 
extremas: 



AT. <vT v 
k X 



0<&T k <(2° -\)T X 



(9a) 



(9b) 



La forma (9a) resulta util para realizar la focalizacion dinamica a partir de una 
distancia minima R mint esto es, para R > Rmin, mientras que la forma (9b) seria 
preferible para focalizar dinamicamente en las regiones mas proximas al array, 
esto es, para R < Rmax. Entre uno y otro caso pueden darse rangos intermedios 
mas adaptados a una aplicacion particular, donde R m i n <R< Rmax. 

Sin embargo, con la Tecnica de Correction Focal Progresiva se prefiere la 
forma (9a) que no limita la maxima distancia de validez. Para las regiones mas 
proximas al array en las que no es posible verificar la desigualdad de la 
izquierda de la Ecuacibn (9a), la tecnica propuesta proporciona una serie de 
altemativas que seran descritas con detalle mas adelante. Asi, pues, llamando 
en la Ecuacion (9a): 



a = 2°-\ 

puede escribirse en forma mas compacta: 

(v-a)T x <AT k <vT x 



(10) 



(11) 



La Ecuaci6n (6) siempre verifica la desigualdad de la derecha de la Ecuaci6n 
(11), de acuerdo con la Ecuacion (8). Para que tambien verifique la de la 
izquierda: 



1 + 



R-x t send 



1 +*\- 2Rx k senQ 



vT 



(12) 



Con: 
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2 + x\ - 2Rx,senQ =(#- x.senQ j 



1 + 



*f cos z 9 



(13) 



la Ecuacion (12) queda como: 



1 + - 



1 



{ i 



i+ 



2 2 
x, cos G 



[p.-x,senof- 



— >v - a 
2 



(14) 



donde, haciendo los cambios de variable: 

Z k = X k ° 0S ^ 

se obtiene: 



1 + 



r = R - x, senQ 
k 



(15) 



\ v k ) 

Despejando rde la ecuacion anterior se obtiene: 

z k (v -2a) 



v _ 

— >v - a 
2 



(16) 



v>2a 



2^av -a 2 

10 que, deshaciendo los cambios de variable de las Ecuaciones (15) conduce a: 



(17) 



R>R Q (x k ,Q) = 



v -2a 



l^av -i 



— X, cosG 
2 h 



XrSenQ 
k 



v>2a 



(18) 



donde |-| representa el valor absoluto y considera el caso mas desfavorable de 
signos relativos de x k y 9. La ecuacidn (18) representa una cota superior para 
la distancia minima Ro(x ks Q) a partir de la cual los intervalos temporales para 
15 muestrear la serial recibida por el elemento k con un angulo de deflexion 9 
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pueden representarse con un numero entero de b bits, con un error inferior a 
medio periodo T x de un reloj maestro. 

Si la Ecuacion (18) se verifica para todos los elementos del array a todas las 
distancias y dngulos de la region a inspeccionar, el retardo inicial para adquirir 
5 la primera muestra, correspondiente al foco Fo, en el canal k viene dado por: 



donde R 0 = Ro(x ki Q) dado por la Ecuacion (18). Obviamente la primera 
adquisicion puede realizarse a una distancia Ri > Ro, en cuyo caso el valor de 
T 0 (x k , Q) se calculara para este nuevo valor. 

10 A partir de este instante, los intervalos de muestreo entre focos presentan un 
error acotado a medio periodo del reloj maestro, pudiendo codificarse para 
cada foco F, de coordenadas (Ri, Q) en el canal k como: 



donde [-]u representa la funcitin de redondeo al entero m£s pr6ximo y ATm (R h 
15 Q) esta dado por el valor de la Ecuacion (6). Sustituyendo en (20) las 
Ecuaciones (10)y (11): 



que indica que los intervalos de muestreo pueden codificarse con b bits, 
produciendo un error inferior a medio periodo de reloj maestro. El valor Q w es el 
20 codigo de correccidn focal que determina el numero de periodos de reloj 
maestro que debe adelantarse la adquisicion de la serial correspondiente a I 
foco / en el canal k. 

La adquisicion de un numero m > 1 de muestras entre focos aumenta la 
eficiencia de memoria al compartir un unico cbdigo de correcci6n focal de b bits 
25 entre varias muestras, con una densidad de b/m bits por muestra. Si, por 
ejemplo, b=1 y m=4, se requieren 0 f 25 bits para focalizar cada muestra 




(19) 



0<2 <2"-l 



(21) 
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adquirida. Una buena eleccion para reducir el error en los intervalos de 
muestreo consiste en distribuir uniformemente el adelanto expresado por el 
c6digo de correccibn focal Q*, entre las m mue stras. Para esto, en una 
realizacion preferida, se calculan los valores: 



*U>q) = 



\<j<m (22) 



donde g^Q* representa, genericamente, el codigo de correccion focal para 
simplificar la notacion y [-]u es la funcion de redondeo. El valor de u(j, q) 
representa el entero mas proximo del numero de periodos del reloj maestro 
entre el primer foco y la muestra ycon j=1, 2, m t consiguiendo una 
1 0 distribucion lo mas uniforme posible del adelanto q entre las m muestras. 

El intervalo temporal entre un foco cualquiera con la muestra que le sigue en el 
orden j-m es u(m, q) = v-g, esto es, la muestra j-m se adquiere exactamente 
en la posicion del siguiente foco, con un adelanto de q ciclos del reloj maestro 
respecto del intervalo nominal v entre focos, de acuerdo con lo expresado por 
1 5 el codigo de correccion focal. 

El intervalo entre muestras, en periodos del reloj maestro, viene dado por: 

An(l, q) = w(l, q) j = 1, por definition 

Aw(y, q) = u(j, q)-uU~U) 2<j<m < 23 > 

En el caso m£s habitual en el que el numero de muestras entre focos es mayor 
o igual al m£ximo adelanto a aplicar entre focos, esto es, q < < m, al 

20 hacer una distribucion uniforme de q ciclos entre m muestras, el adelanto gj 
asignado a la muestra j solo puede tomar los valores 0 6 1 respecto al intervalo 
nominal v/m que se da para g=0. Este hecho facilita la realizacion de la 
distribucion uniforme del adelanto q entre un numero de muestras m >q 
mediante una tabla en la que se codifica con 1 bit el valor g(j,q), previamente 

25 calculado segun: 

gU>q)=v/m-AuU,q) \<j<m y 0<q<2 b -l (24) 
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Para clarificar el metodo, considerese un ejemplo en el que m=3, v=12=4 m, 
b=2. Los posibles valores de q son 0, 1 , 2 y 3 por estar codificado con 2 bits, y 
v/m=4. La Tabla I muestra los resultados que se obtienen al aplicar, con estos 
valores, las Ecuaciones (22), (23) y (24). Se observa que el adelanto qse 
5 distribuye uniformemente entre las muestras dentro de la resolucion de 1 
periodo de reloj maestro para todos los codigos de correction focal. La 
densidad de la codificacion en este ejemplo es de 2/3 = 0.67 bits/muestra. 

Tabla I. 

10 Ejemplo de distribution uniforme de q adelantos entre m muestras 





v-q 


u(1.q) 


u(2,q) 


u(3,q) 


Au(1,q) 


Au(2,q) 


Au(3,q) 


9(1 .Q) 


9(2,q) 


9(3,q) 


q=0 


12 


4 


8 


12 


4 


4 


4 


0 


0 


0 


q=1 


11 


4 


7 


11 


4 


3 


4 


0 


1 


0 


q=2 


10 


3 


7 


10 


3 


4 


3 


1 


0 


1 


q=3 


9 


3 


6 


9 


3 


3 


3 


1 


1 


1 



Para cada valor de q hay m valores posibles gQ,q)> cada uno de 1 bit. Por 
consiguiente, todos los adelantos individuates pueden codificarse en una 
memoria de m ^ palabras de 1 bit (en el ejemplo anterior, basta una memoria 
15 de 12x1 bit). 

Para valores R < Ro(x k , Q) el error absolute del intervalo entre instantes de 
muestreo puede superar el valor 7>/2. En este caso, las tecnicas CFP y CFPV 
pueden operar de varias formas, destacando las siguientes estrategias: 

1. Activar una funcion de apertura dinamica en la que el elemento /cno 
20 intervene hasta que se reciban senales desde una distancia igual o superior 
a Rofok, Una de las posibles realizaciones, produce para este elemento 
un numero de muestras nulas M z (k) equivalente a las que se producen a la 
salida del conformador entre el comienzo de la adquisicion T ini y el instante 
To(x ki d) dado por la Ecuacion (19): 
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\(x k ,e)-T. 



M z ik) = 



mi 



(25) 



donde [■]? representa la funcion de redondeo por exceso a un numero 
entero. 

Debe destacarse que la relation D/max(R 0 ), donde D es el tamano de la 
5 apertura, determina el mmimo valor de F# que puede aplicarse con las 
t6cnicas propuestas. En un ejemplo prSctico se muestra que F# m in =1, 0'6 y 
0'4 para valores de b= 1, 2 y 3, respectivamente, muy competitivos con los 
proporcionados por otras tecnicas. 

2. Anticipar la adquisicion de senal en el elemento k a una distancia Ri < 
10 Ro(*k, elegida de forma tal que se adquiera una muestra en R 0 (x k , Q). Los 

instantes de muestreo anteriores a este tendran un cierto error superior al 
acotado, pero su efecto en la imagen puede ser limitado. Para ello, el 
instante de adquisicion de la primera muestra dado por la Ecuacion (19) 
queda modificado a: 

15 r Q M (x k ,6 ) = T Q {x k , 9 ) - M z (k)T sl (26) 

donde se expresa que, para el elemento k t se adquieren M z (k) muestras 
anteriores a la reception con error acotado y Tsi es el periodo de muestreo 
utilizado, constante en esta region. El valor de T S i puede ser igual al 
periodo de muestreo de salida del conformador 7s o, preferiblemente, un 
20 valor inferior para reducir en lo posible los errores en el instante de 
muestreo. 

3. Operar de forma mixta, donde un primer rango de distancias [R 2 , Ri) se 
cubre con apertura dinamica, otro [Ri, Ro) con un error limitado como en el 
punto anterior y, finalmente, un tercero [Ro, <*>) en el que los errores de 

25 muestreo estan acotados dentro de medio periodo de reloj maestro 
aplicando la correction focal progresiva. 

En una realization preferida que, ademas resulta de aplicacion general cuando 
la geometna del array es diferente, el medio no es homogeneo (por ejemplo, 
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cambios de velocidades de propagation) o los focos se distribuyen a intervalos 
regulares pero siguiendo una trayectoria diferente a la considerada, se utiliza 
un algoritmo de cSlculo basado en los siguientes pasos: 

ALGORITMO 1: 

5 PAS0 1 : Se calculan los tiempos de propagacitin T ki de la senal 
correspondiente a cada foco F/, /= 0, 1, 2, hasta cada elemento k 
del array. El valor de Tu se obtiene en cada caso al considerar la 
geometria del sistema y la(s) velocidad(es) de propagation. 

PASO 2: Para cada elemento k, se calculan los intervalos U ki entre instantes 
10 de llegada de la senal correspondiente a dos focos consecutivos, 

expresados en periodos del reloj maestro, como: 



T - T 

'ki it, i — 1 



T 



i>\ (27) 

n 



PASO 3: El valor v o intervalo nominal en periodos de reloj maestro entre 
focos se elige de los resultados de la ecuacion anterior como: 

15 v>max(C/ fa .) V*,i (28) 

PASO 4: Si el niimero de bits disponibles para codificar las correcciones 
focales es b, determinar para cada canal k el indice i=h k a partir del 
cual: 

v-U H <2 b -\ i>h k (29) 

20 PASOS: Para los focos F/con /< h kl el numero de muestras M z (k)> para 
operar con apertura din£mica o adquisicibn anticipada en el canal k t 

es: 

M z (k)=(h k -\)m (30) 

PASO 6: Los codigos de correction focal para cada elemento k 
25 correspondientes al foco / > h k se calculan como: 

Q u =v-U M i>h k (31) 
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que expresan el numeric de ciclos de reloj maestro que debe 
adelantarse la adquisicion de la senal correspondiente al foco / en el 
canal k. 

En otra realizaci6n preferida, durante la adquisicion de senal, la distancia entre 
dos focos consecutivos F h1 y F\ se varia en momentos seleccionados 
incrementando el numero de muestras que se adquieren entre ellos. Esta 
variante se denomina Tecnica de Correction Focal Progresiva Variable o 
CFPV, en la que el numero de muestras adquiridas entre dos focos 
consecutivos pasa de ser una constante m a una funcion del numero de foco 
actual, esto es M=m(i). 

La Ecuacion 18 indica que el valor de R 0 crece con el de v. Por su parte, el 
valor mmimo de v esta determinado por la Ecuacion 1 cuando se adquiere una 
muestra por foco, con un intervalo temporal entre muestras Ts = sT*. Iniciando, 
pues la adquisici6n con M=m(1)=1 a partir de una distancia R 0 (1) determinada 
por el valor v = s a la que, para todos los elementos: 



s-2a 



Was -a 2 



JCj, COS0 



x^senQ 



(32) 



donde a = 2 b - 1 . La adquisicidn continua manteniendo M=m(1)=1 hasta 
alcanzar cierto foco a una distancia R=R 0 (2) a partir de la cual se verifica, para 
todos los elementos: 



20 



* 0 (2)> 



2s -2a 



> V 2as - a 



x, cos9 
2 k 



x, send 
k 



(33) 



25 



donde comienza un tramo en el que pueden adquirirse 2 muestras por foco, 
esto es, M=m(2)=2 t manteniendo el error de muestreo acotado. Cierto tiempo 
despu6s, se alcanzarS otro foco en R=Ro(3) a partir del cual puede operarse 
con M=m(3)=3 y as! sucesivamente, de forma que, genericamente, a partir de 
un foco situado en R=R 0 (m) a partir del cual pueden adquirirse m muestras, 
verificandose: 
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ms — 2a 



jc, cos9 + x 1 send Vfc 



(34) 



2l 



ams 




En definitiva, con la Tecnica de Correccion Focal Progresiva Variable el numero 
de muestras que se adquieren entre focos se incrementa con el tiempo, de 
forma que los codigos de correccion focal puedan seguir expresandose con b 
5 bits. Debe destacarse que la condition (34) debe ser verificada por todos los 
elementos para cada valor de m. 

En una realization preferida, se asocia a cada foco F/ un campo de 1 bit, 
denominado codigo de alejamiento focal J/. Cuando JfO se mantiene el actual 
valor m(i)=m(i-1) y, cuando J,=1 se incrementa: m(i)=m(i-1)+1. El algoritmo 2 
10 proporciona un posible metodo para el calculo de los valores de Jf. 

ALGORITMO 2: 

PASO 1 : En funcion de la distancia minima establecida por la aplicacion, R mirh 
determinar el mayor valor m=m 0 para el cual se verifica la Ecuacion 



PASO 2: Inicializacion: m=m 0 . Variables auxiliares: f=1 (numero actual de 
foco), /=7 (numero actual de muestra). 

PASO 3: Repetir: 

a) Calcular R 0 (m) dado por la Ecuacion (34) y el valor AR = cTs/2, 



15 



(34). 



20 



siendo c= velocidad de propagacion, T s = periodo de 



muestreo. 



b) Si R min + (i+m) AR > (m+1)R 0 



{ Jf- 1; m=m+1;} 



en otro caso {J r = 0;} 



25 



c) tt+1; 



d) i-hm; 



hasta que / > A/s, siendo Ns = numero de muestras a adquirir. 
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El numero de codigos de alejamiento focal iguala al de focos para cuaiquier 
adquisicion, por lo que puede realizarse facilmente asociando a cada codigo de 
correccibn focal Q M el correspondiente cbdigo de alejamiento focal J w = J t en 
cada canal k. 

5 Este esquema resulta adecuado ya que, en general, las senales 
correspondientes a cada foco llegan a los elementos en momentos diferentes y 
su realizacion es sencilla. De este modo, la memoria de correccibn focal tiene 
una anchura de palabra de b+1 bits, repartidos en dos campos: b bits para el 
codigo Qwy 1 bit para el J h 

10 Una alternativa, que proporciona una mayor eficiencia, almacena en una 
memoria separada los codigos de alejamiento focal J/, cuyo contenido es 
compartido por un conjunto de canales. Para resolver el problema de no 
simultaneidad en la adquisicion de los focos, cada canal dispone una pequena 
cola donde se introducen los codigos J, obtenidos autom£ticamente de la 

15 memoria que los almacena. Cada canal utiliza el valor de J/ en cabeza de la 
cola cada vez que alcanza un nuevo foco. Asf, una unica memoria central 
almacena un codigo J, por foco y unos pequenos FIFOs distribuidos facilitan la 
asignacibn de estos cddigos comunes a cada canal cuando lo necesita, como 
se expone con mas detalle en la description del sistema. 

20 Por otra parte, puesto que la distancia entre focos puede modificarse 
dinamicamente, en la mayoria de las aplicaciones practicas basta una 
codificacion de la correccion focal con b=1 bit, lo que proporciona una alta 
eficiencia en el uso de memoria y una simplification de los circuitos de 
correccibn. En particular, es inmediato realizar la distribucibn uniforme de 0 6 1 

25 periodo de reloj maestro entre m muestras. 

Con lo expuesto, la Tecnica de Correccion Focal Progresiva, con su opcion 
Variable, proporciona un metodo que realiza la conformacion dinamica de 
haces en recepcion con una elevada resolucion temporal dentro de un gran 
rango de distancias. Adicionalmente, carece de redundancias en la adquisicion, 
30 no requiere procesos de interpolacion de muestras, permite mantener la tasa 
de muestreo mas baja compatible con el criterio de Nyquist y, por consiguiente, 
tiene un bajo consumo energetico. Su realizacion es compacta, debido tanto a 
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un balance adecuado entre el calculo en tiempo real de las conrecciones y el 
contenido de unas memorias precalculado en tiempo no real, como a que la 
codificaci6n de esta informaci6n es muy eficiente. 

Todo esto pemnite configurar un sistema vers£til de conformacion de haces en 
5 recepcion con focalizacion dinamica, deflexion, apertura dinamica y 
apodizacion, cuyos detalles constructivos se exponen a continuation. 

Description del sistema. 

El sistema objeto de la presente invention se basa en la tecnica de Correction 
Focal Progresiva descrita anteriormente, incluyendo o no la option Variable, y 
10 realiza la focalizaci6n dinamica con deflexion de las senales recibidas por un 
conjunto de transductores o array, incluyendo funciones de apodizacion y de 
apertura dinamica. El sistema se describira con suficiente detalle para que 
pueda ser entendido y reproducido por personal con conocimientos de 
electronica digital. 

15 Con referenda a la Figura 2, un sistema de conformacion de haces en 
recepcion consta, en el caso mas general, de un array (10) compuesto de N 
transductores elementales (11) o, simplemente, elementos, que reciben 
senales del exterior, numerados 1, 2, .../c, N. Frecuentemente, estos 
elementos son, a su vez, emisores de senal para conformar un haz en emision, 

20 pasando a ser receptores una vez que esta finaliza. En la figura se han omitido 
los dispositivos necesarios para conformar el haz en emision. 

Las senales recibidas por los diferentes elementos son acondicionadas y 
amplificadas por los dispositivos AAS (12), que entregan las N senales 
analogicas ai, a* a*> aw a otros tantos conversores analogico-digitales 

25 A/D (13). Los dispositivos AAS suelen incorporar, ademSs de la funci6n de 
amplificaci6n, otras que incluyen la compensaci6n de la atenuaci6n con la 
distancia mediante un control ganancia-tiempo, el filtrado y, en ocasiones, el 
desdoblamiento en las componentes en fase y cuadratura. En este ultimo caso, 
cada senal analogica a* esta formado por el par (/*, q k ) que representan, 

30 respectivamente, a las componentes en fase y en cuadratura. 
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Una vez acondicionadas, las senates analdgicas son digitalizadas por los 
conversores A/D. En el caso de senales en banda base con pares (/*, q k ) 
pueden utilizarse dos conversores A/D por cada senal o multiplexarse un unico 
conversor A/D entre dos senates, las muestras pares para la componente i k y 
5 las impares para la q k o viceversa. 

Otra alternativa es la utilization de conversores AZ, con los que las senales 
analogicas a* se representan mediante una secuencia continua de bits. En este 
caso la tecnica de Correccion Focal Progresiva y sus variantes resulta 
especialrnente util al ajustar los instantes de muestreo a los de llegada de la 
10 senal a cada elemento con la resolution de un reloj maestro cuyo periodo se 
adapta, ciclo a ciclo, al requerido por tales conversores AL para obtener una 
determinada relation senal/ruido. 

Las tecnicas mencionadas son conocidas por los expertos en el tema y de 
aplicacidn general. Las diferencias aparecen principalmente en el dispositivo 
15 CONF (14) o conformador de haz, cuya funcion es la de combinar las N 
secuencias numericas 9i, e 2l e k , ew proporcionadas por los conversores 
A/D para crear una unica secuencia de salida r que representa, en forma 
digital, a la serial recibida por el array focalizada, idealmente, en todos los 
puntos de una direction determinada. 

20 Cada muestra entregada por los conversores A/D es una palabra de w bits 
siendo, tipicamente, 6 < w < 12 para el caso de conversores de muestreo 
instantaneo y 1 < w < 2 para los de tipo AZ, aunque el sistema puede operar 
con otras anchuras de palabra. 

Es, precisamente, la forma en que el dispositivo CONF realiza la funci6n de 
25 combinar coherentemente las N entradas e k para generar la salida r lo que 
constituye el aspecto diferencial de la presente invencion. 

El dispositivo CONF (14) que es objeto de la presente patente, basa su 
funcionamiento en la tecnica de Correccion Focal Progresiva que constituye 
uno de sus aspectos claves, aportando diversas ventajas frente a otros tipos de 
30 realizaciones. 
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Asf, segun otros metodos resenados en la description que se hizo del estado 
del arte, cada muestra de la secuencia numerica suministrada por k-esimo 
conversor A/D es sometida a un proceso de retardo como indica la Figura 3. 
Cada elemento se asocia a un canal de procesamiento individual (21). El 
5 contenido de una memoria de focalizacion dinamica (25) es utilizado por un 
dispositivo de control (24) para generar senales que producen un retardo 
grueso (22) y otro fino (23). El retardo grueso tiene, tipicamente, la resolucibn 
de un periodo del reloj de muestreo y el retardo fino una fraction del mismo. 
Estos retardos se obtienen con diversos metodos que han sido descritos en la 
1 0 literatura especializada y son bien conocidos. 

El contenido de la memoria de focalizacion dinamica (25) es tal que, para cada 
muestra y para cada canal, los retardos introducidos alinean en el tiempo las 
salidas f 1t fz f/v que corresponden a los focos sucesivos F 7 , F 2 , F N . Estas 
senales son sumadas en (20) para obtener la secuencia r focalizada 
15 dinamicamente. 

La mayor ventaja de estas tecnicas reside en que solo es necesario un reloj de 
muestreo, comun a todos los canales. Pero adolece de diversos inconvenientes 
derivados de la necesidad de realizar procesos de interpolacibn para obtener 
una suficiente resolution temporal en los retardos introducidos a las senales, 
20 con el consiguiente aumento de complejidad electronica y la posibilidad de 
introducir errores. Ademas, requieren una considerable cantidad de memoria 
para representar los retardos grueso y fino. 

El sistema objeto de la presente invencibn opera de una forma diferente, 
muestreando directamente las senales en el instante en que llegan a cada uno 

25 de los elementos del array, procedentes de cada uno de los focos situados a lo 
largo de una direction determinada. Al evitar la introduction de retardos a las 
senales, se hacen innecesarios los procesos de interpolation, eludiendo 
posibles errores y eliminando los circuitos electrdnicos correspondientes. 
Asimismo, solo se adquieren las muestras necesarias para realizar la 

30 composition coherente, eliminando redundancias. 

En principio, el coste del metodo objeto de esta patente es que cada conversor 
A/D debe operar con un reloj de muestreo diferente, de frecuencia no uniforme. 
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Como se demuestra con la descripcion que sigue, la realization de esta 
caracteristica es sencilla a partir de la metodologia enunciada como tecnica de 
Correcci6n Focal Progresiva, con la opci6n Variable, pudiendo integrarse 
multiples canales en un unico dispositivo tipo FPGA estandar y comercialmente 
5 disponible. 

Con referenda de nuevo a la Figura 2, el metodo propuesto requiere que el 
dispositivo de conformation de haz CONF (14), proporcione un reloj de 
muestreo independiente y de frecuencia no uniforme CKi, CK 2 , CK k , CK N 
a cada conversor A/D (13), de forma que se adquieran las senales analogicas 

10 ait 32, Ski 9n proporcionadas por los acondicionadores de senal AAS 
(12) en los instantes que corresponden a los de su llegada desde cada foco a 
cada elemento (11) del array (10). Las muestras obtenidas por los N 
conversores A/D (13) son procesadas para obtener su suma coherente, es 
decir, un valor de salida r que es la suma de las N muestras correspondientes a 

1 5 la serial recibida de un mismo foco por los N elementos del array. 

En la presente invention, esta metodologia se aplica en forma modular, 
asociando un canal de procesamiento a cada elemento del array. La Figura 4 
muestra la estructura de uno de los canales (30) que procesa la senal recibida 
por el elemento k (11) del array. La estructura de todos los canales es identica 
20 a la mostrada. 

Como en otros casos, la senal es amplificada y acondicionada por el dispositivo 
AAS (12) por metodos conocidos, que incluyen dispositivos electronicos para 
amplificar, fiitrar y manipular de forma anal6gica la senal. A la salida del AAS se 
obtiene la senal a k que es digitalizada por el conversor A/D (13). Las muestras 

25 son adquiridas en los instantes definidos por un generador de reloj (70) que 
produce una secuencia de pulsos CK k no uniformemente espaciados, cuya 
descripcion detallada se realiza posteriormente. El generador de reloj es 
habilitado por la senal HFD generada en el dispositivo HFOC (80) que marca 
los instantes en los que el canal comienza y termina la adquisicion de senal con 

30 focalizacion dinamica. 

La funcion del generador de reloj es la de garantizar que los instantes de 
muestreo determinados por CK k corresponden a los de la llegada de la senal 
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desde cada uno de los focos al elemento k 9 y opera segun los principios 
descritos anteriormente como Tecnica de Correccion Focal Progresiva, con su 
posible opcibn Variable. 

Con esto, la secuencia de muestras e k obtenida a la salida del conversor A/D 
5 corresponde, precisamente, a los valores ordenados de la serial recibida por el 
elemento k procedente de cada uno de los focos. Estos valores son 
procesados por un bloque de apodizacion y apertura dinamica APD (60), que 
permiten mejorar las caracteristicas del haz y que produce la secuencia f k . Al 
existir una correspondence 1:1 de las muestras e k a las f kt esta secuencia 
10 tambien esta ordenada con referenda a los instantes de llegada de la serial al 
elemento k desde cada foco. Los sucesivos valores de f k se almacenan 
transitoriamente en una cola o memoria FIFO que absorbe la diferencia en los 
instantes en que se producen muestras correspondientes a un mismo foco en 
diferentes canales. 

15 Elemento clave es, pues, el generador de reloj de muestreo (70), que produce 
la serial de digitalizacion CK k del k-esimo conversor A/D, de forma que las 
muestras obtenidas correspondan a los valores de la serial recibida por el 
elemento k desde cada uno de los focos. La Figura 5 muestra el diagrama de 
principio de una de las realizaciones preferidas del generador del reloj CK k , 

20 cuyo funcionamiento se basa en el metodo descrito anteriormente como 
Tecnica de Correccion Focal Progresiva, para el caso en que las correcciones 
focales se codifican con 6=1 bit. 

El circuito mostrado en la Figura 5 contiene s6lo los elementos bSsicos 
omitiendose otros dispositivos auxiliares, que son necesarios para el correcto 

25 funcionamiento, en favor de una mayor claridad de exposition. Tiene un 
registro REG-A (71a) en el que se carga un valor relacionado con el parametro 
v o intervalo nominal entre focos en periodos del reloj maestro. Su salida V 
consta de multiples h'neas que determinan el numero de etapas que componen 
el registro de desplazamiento de longitud programable SHR-A (72). Este 

30 dispositivo obtiene una salida p que es una versi6n retrasada una cantidad de 
V periodos de reloj c de su entrada CK k . tzste y otros dispositivos en el 
esquema mostrado operan sincronamente con el reloj maestro c de periodo T x . 
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La serial de entrada HFD habilita la operation del registro de desplazamiento 
(72) y del flip-flop (77), de forma que estos dispositivos mantienen su estado 
actual independientemente del valor a su entrada cuando HFD=0. Cuando 
HFD=1, se habilita la focalizacion dinamica, interrumpiendose cuando HFD 
5 retorna a cero. Esta serial se genera en el dispositivo HFOC (80), Fig. 4. 

La senal p es, a su vez, retrasada 1 periodo de reloj maestro mediante el flip- 
flop FFA (73), produciendo la senal pr. El multiplexor MUX (74) proporciona una 
salida s que es, bien la senal p producida por el registro de desplazamiento o 
su version retrasada pr, en funcion del estado de la Ifnea de selection Q: 

10 s=p si Q=0; s=pr siQ=1 (35) 

El estado del selector Q lo determina el contenido de una memoria MEM (75) 
en la direction F proporcionada por el contador CNT (76). En el caso actual, Q 
es una simple linea (1 bit). 

La memoria se carga previamente, mediante circuitos no mostrados, con los 
15 codigos de correction focal Q*„ donde k es el numero del canal actual y / es el 
numero de orden del foco, con /> 1. Los circuitos para la adquisicion de senal 
en el primer foco (/=0), que se describen posteriormente, determinan la 
activacidn de la senal HFD. 

La salida del multiplexor (74) se registra en el flip-flop FFB (77), cuya salida 
20 CK k es el reloj de muestreo del conversor A/D del canal actual k (13), Figura 4, 
acciona al contador CNT (76) y es la entrada del registro de desplazamiento 
SHR-A(72). 

Para su funcionamiento, el registro de desplazamiento SHR-A (72) se inicializa 
a '0\ el flip-flop FFA (73) a '0\ el flip-flop FFB (77) a T y el contador (76) con la 
25 direction de memoria donde se almacena el codigo de la correction focal 
correspondiente al foco F 1 (i=1). Los circuitos auxiliares para establecer estos 
estados iniciales se omiten en la figura. Asi, los valores iniciales de las 
diferentes senales son: CK k =1, p=0, pr=0 ys=0 (independientemente del 
estado de Q). 

30 A partir de este estado initial, en el ciclo de reloj maestro c que sigue a la 
activation de la senal HFD=1, se pone CKk=0 al encontrarse en este estado la 
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entrada s del flip-flop FFB (77) y el reloj habilitado por su entrada ce. Durante 
los V primeros periodos del reloj c se mantienen las senales p=0 y pr=0. En el 
ciclo siguiente se pone p=1, reflejando que el valor inicial CK k =1 se ha 
propagado por el registro de desplazamiento (72), pero todavia se mantiene 
5 pr=0. En este mismo ciclo, si Q=0, es s=pr=0, pero si Q=1 , es s=p=1. Por 
consiguiente, en el siguiente ciclo, la salida CK k =1 solo si Q=1\ tambien en 
este ciclo se pone pr=f por lo que, si Q=0 se pondra la salida CK k =1 en el ciclo 
siguiente. Una vez que se hace CK k =1 el proceso se repite. 

Por consiguiente, el numero de ciclos n que transcurren desde el anterior valor 
10 CK k =1 y el proximo lo determina el valor Q de la siguiente forma: 

Si Q=0: n=V+1 

S\Q=1: n=V+2 (36) 

Acorde con la Ecuacion (20), cuando Q=0 el numero de ciclos entre focos es n 
= v y, cuando Q=1, debe ser n=v-1. Con esto: 

15 V=v-2 (37) 

que indica el valor inicial a cargar en el registro REG-A (71a) en funci6n de la 
distancia nominal entre focos v expresada en periodos del reloj maestro. 

Cada vez que se produce un reloj de muestreo (CK k =1 ), se incrementa la 
direcci6n F proporcionada por el contador CNT (76), con lo que la memoria 
20 MEM (75) coloca en su salida el c6digo correspondiente a la correcci6n focal a 
aplicar en el siguiente foco. 

Se observa que el conjunto FFA (73)-MUX (74)-FFB (77) configuran un registro 
de desplazamiento de longitud variable segun el codigo Q, que puede utilizarse 
como base para las realizaciones descritas a continuacion, aunque hay 

25 altemativas m£s eficientes. Asf, la Figura 6 muestra un esquema de otra 
realizacion preferida que no limita el numero de bits con el que se representan 
los codigos de correccion focal Q. En este caso el conjunto FFA-MUX-FFB 
anterior se ha sustituido por un registro de desplazamiento SHRB (78) cuya 
longitud programable la determina el complemento del codigo de correccibn 

30 focal Q obtenido mediante los inversores (75a) que producen las salidas A/Q. 
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En este caso el intervalo entre los instantes de muestreo en dos focos 
consecutivos en el canal actual se determina.por el valor multibit Qentregado 
por la memoria (75). Aqui el registro de desplazamiento SHRB (78) se inicializa 
aT. 

5 El funcionamiento es, por lo demas, analogo al caso anterior con la salvedad 
de que el numero de ciclos esta determinado por los valores V y NQ segun: 



Teniendo en cuenta que A/Q=2*-7-Q, el valor inicial a cargar en el registro REG- 
A (71a) en funcion del intervalo nominal entre focos en periodos del reloj 
10 maestro es: 



En una realization practica, los inversores (75a) son innecesarios si en la 
memoria de codigos de correccion focal (75), se almacenan los valores 
complementarios NQ en lugar de los de Q, manteniendose en la description 
15 por coherencia en la notation. Asimismo, en las realizaciones practicas de los 
circuitos mostrados en las Figuras 5 y 6 se incluye, ademds, I6gica auxiliar que 
evita la situation en la que todas las salidas estan en estado '0' de la que, 
evidentemente, no saldrian los circuitos. 

La logica descrita permite focalizar cada una de las muestras adquiridas. 

20 Cuando se prefiera realizar una correccion focal periodica, compartiendo un 
mismo c6digo de correccibn focal Q cada m muestras, los circuitos anteriores 
se modifican como indica la Figura 7. El valor de m se carga en un registro 
RCM (71b) y el valor almacenado M es utilizado por una logica de control 
CTRL (90) junto con el valor Q del codigo de correccion focal entregado por la 

25 memoria (75) para generar las senales ng y cef. La primera controla la longitud 
del registro de desplazamiento SHRB (78) haciendola igual a ng+f , mientras 
que cef habilita al contador CNT (76) que proporciona la direccion F a la 
memoria de correccion focal (75). 

El registro de desplazamiento SHRA (72) es de longitud programable en 
30 funcibn del valor V programado en el registro REG-A (71 a), con un numero de 



n=V+NQ 



(38) 



V=v-2 b +l 



(39) 
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etapas igual a V+f. El registro de desplazamiento SHRB (78) es tambien de 
longitud programable con un numero de etapas igual a ng+f . 

El numero de periodos de reloj maestro que median entre dos ciclos 
consecutivos CK k es igual al numero de etapas interpuestas entre la entrada de 
5 SHRA (72) y la salida de SHRB (78) que, acorde con lo anterior, es: 

n = V + ng + 2 (40) 

Programando en el registro REG-A (71a) un valor 

K=i)/m-3 (41) 

donde v/m representa el intervalo nominal entre muestras en periodos del reloj 
10 maestro, la Ecuacion (40) queda como: 

n =v /m + ng-l (42) 

La logica de control CTRL (90), mostrada en detalle en la Figura 7, controla el 
intervalo n mediante la serial ng, distribuyendo uniformemente el adelanto de 
periodos de reloj maestro Q entre las m muestras que se adquieren entre dos 
1 5 focos consecutivos. 

Para ello, el contenido CM de un contador m6dulo m, CMOD (91), que se 
actualiza con cada reloj CK k por su entrada de habilitacion ce, proporciona 
parte de la direccion DM de la memoria MDM (92). La otra parte de la direccion 
es el codigo de correction focal Q. La direccion proporcionada a la memoria 
20 MDM es: 

DM(j,QHj-m b +Q (43) 

donde /= CM+1 indica el numero de orden de la muestra actual. Programando, 
en tiempo no real, el contenido de las direcciones DM(j,Q) de la memoria MDM 
con los valores g(j,Q) que resultan de aplicar las Ecuaciones (22) a (24), la 
25 salida g representa el adelanto individual a aplicar a cada una de las muestras 
durante la adquisicion. El inversor (93) obtiene la serial ng, complementaria de 
la g, de modo que, conforme a la Ecuacion (42): 

Si g=0 => ng=7 n = v/m 
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Si g=1 => ng=0 => n = v/m - 1 (44) 

con lo que se produce un adelanto de 1 ciclo de reloj cuando g(j,Q) =1 y se 
mantiene el intervalo nominal v/m cuando g(j,Q) = 0, tal y como qued6 descrito 
en el m6todo de Correcci6n Focal Progresiva. 

5 Por su parte, cada vez que el contador CMOD (91) rebosa (esto es, conmuta 
del valor m-ta 0), se produce un pulso cef que habilita el increment© del 
contador de focos CNT (76), proporcionando una nueva direccion F+7 a la 
memoria MEM (75). Asi, el proceso se repite con un nuevo valor de Q 
suministrado por la memoria MEM (75) con el contenido de la nueva direccitin. 

10 Deben hacerse algunas observaciones adicionales. La primera es que, 
evidentemente, los inversores (93) se eliminan si en la memoria se programan 
los valores complementarios de los g(j,Q) dados por la Ecuacion (24). Otra es 
que la memoria MDM (92) es de pequeno tamano; por ejemplo, una memoria 
de 256x1 bits puede servir para valores m< 64 con b=2, o para m< 32 con b=3. 

15 A partir de este esquema se deriva la variante denominada como tecnica de 
Correction Focal Progresiva Variable. Para ello, el registro RCM (71b) se 
sustituye por un contador, de forma que el numero de muestras m entre focos 
puede irse incrementando en el tiempo a partir de un valor mmimo inicial m 1 
que, ti picamente, puede ser igual a 1 . 

20 Tal y como fue descrito, segun la Ecuacion (18), la distancia minima R 0 a partir 
de la cual el error absoluto del instante de muestreo se mantiene acotado en 
±Tx/2 crece con el valor del intervalo entre focos v, expresado en periodos T x 
del reloj maestro. Al reducir la distancia entre focos en la region mas proxima al 
array mediante la programacion de un valor bajo para m* f disminuye el valor de 

25 R 0 . Los instantes en que se incrementa el contador RCM (71b) vienen 
determinados por los instantes en que las senales proceden de focos situados 
a suficiente distancia como para que pueda adquirirse una muestra mas entre 
focos en todos los canales. El numero de orden de los focos en que esto ocurre 
se calcula en tiempo no real y se codifican en una tabla comun a todos los 

30 canales, a un subconjunto de canales o a un canal individual, estableciendo un 
balance entre eficiencia de uso de memoria y modularidad. 
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En una realization preferida y, en tiempo de adquisicion, los focos en los que 
se incrementa el contador RCM (71b) se determinan por el valor de un codigo 
de alejamiento focal J, expresado con 1 bit, con / = f, 2, ... asociado a cada 
foco: el valor actual de m contenido en el contador RCM (71b) se incrementa si 
5 Ji = 1 , manteniendo su valor actual si J/ = 0 hasta el siguiente foco, donde se 
repite este proceso. El incremento de RCM solo se habilita, por consiguiente, 
cuando se produce el reloj de muestreo para adquirir la serial correspondiente 
a un foco. 

En una realization particular, la memoria de correction focal (75) proporciona 
10 dos campos: uno, J fe = Ji de un bit, indica si se incrementa el contador RCM 
(71b) en este foco y, el otro Q w de b bits, el adelanto a distribuir entre las m 
muestras actuates. En la figura se muestra la senal J=Jj con una li'nea de 
trazos. En este caso, preferentemente, se utilizan c6digos de correccibn focal 
Q w con b-1 bit, puesto que puede controlarse de forma separada la distancia 
15 entre focos controlando el numero de muestras m intercaladas. En tal caso, la 
distribution de los adelantos expresados por Q w son de 0 6 1 ciclo de reloj 
maestro, realiz£ndose en la muestra cuyo orden es la mitad del contenido de 
RCM (71b), aunque son posibles otras alternativas mas simples con poca 
incidencia en la calidad de la imagen. 

20 Sin embargo, para mejorar la eficiencia, en otra realization altemativa que 
tambien se muestra en la Figura 7, los codigos de alejamiento focal Ji 
proceden de una memoria separada (MJ (59), Figura 10) cuya salida es 
compartida por multiples canales por la entrada J a una cola FIFO local (79). 
En tiempo de adquisicion, se realiza una lectura de un nuevo codigo de 

25 alejamiento focal J, en la memoria global MJ (59) de la Figura 10, cuando 
alguno de los canales a los que sirve detecta que su FIFO local (79) este vacio. 
El codigo J, suministrado se escribe simultaneamente en todos los FIFOs 
locales (79) para ser utilizados por cada canal cuando adquieren serial de un 
nuevo foco, esto es, se activa su senal eel En este momento, el valor en 

30 cabeza de la cola Jc se utiliza para incrementar o no el valor de m contenido en 
el contador RCM (71b). 
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Debe destacarse que las colas FIFO locales (79) son de pequeno tamano, 
tipicamente inferiores a 16x1 bit porque solo ha de almacenar un numero de 
c6digos de 1 bit equivalente a la maxima diferencia de focos adquiridos en 
diferentes canales dentro de un mismo modulo. De este modo, con esta 
5 altemativa se consumen menos recursos que con la anterior en la que los 
codigos de alejamiento focal J, se almacenaban en la memoria de correccion 
focal (75) de cada uno de los canales. 

Con la Tecnica de Correccion Focal Progresiva y Variable la utilizacion de 
codigos de correccion focal se hace mas eficiente, distribuyendose entre un 
10 numero bajo de muestras al principio de la adquisicion (en el limite, 1 bit por 
muestra) y entre un numero mayor de muestras conforme la distancia de los 
focos al array aumenta (una fraction decreciente de un bit por muestra). 



Por otra parte, es preciso determinar con precision el instante en el que se 
15 adquieren las muestras que corresponden al primer foco, de coordenadas (R 0 , 
Q) t a partir del cual opera la correccibn focal segun la metodologfa descrita. 
Este instante este determinado por la Ecuacibn (19) para la geometria de la 
Figura 1 . Debe hacerse la observation que el foco inicial puede situarse a una 
distancia mayor en funcion de las necesidades de la aplicacion. Tambien se 
20 tendra en cuenta la utilizacion o no del un numero variable de muestras entre 
focos para determinar este valor inicial. 

Esta funci6n se realiza en el bloque HFOC (80) de la Figura 4, expuesto con 
mas detalle en la Figura 8. Una realization utiliza, simplemente, un contador 
accionado por el reloj maestro que se carga inicialmente con el valor Na(I<): 



25 N A (k) = 



T 0 (k) 



T x Jti 



(45) 



donde T 0 (k) representa el intervalo entre el origen de tiempos hasta que la 
serial correspondiente al foco Fo Hega al elemento /c, determinado por la 
geometria y las velocidades de propagacion, o como expresa la Ecuacion (19) 
para la configuraci6n de la Figura 1. 
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Las senales que corresponden al primer foco F 0 se adquieren, en el canal k, 
cuando rebosa este contador, habilitando, a partir de este momento, la 
actuaci6n de los circuitos de correcci6n focal progresiva. El error absoluto del 
instante de adquisicion de este primer foco es inferior a medio periodo de reloj 
5 maestro en todos los canales. Sin embargo, para mantener esta resolution 
temporal, este contador debe ser accionado a la frecuencia del reloj maestro, 
normalmente elevada y muy superior a la de muestreo. 

Para evitar este inconveniente, una realization preferida, descompone el 
contador de ciclos N A en dos partes, como se muestra en la Figura 8, de forma 
1 0 que (omitiendo el indice de canal k): 

N A =K-N A1 +N A2 + 1 (46) 

con K } N A i y Na2 enteros. El valor N A i se carga en el contador CNT-A (82) 
mediante la serial Ida. Cuando se habilita su funcionamiento mediante la 
activacion.de la serial cea, se acciona el contador con un reloj cb cu ya 
15 frecuencia es K veces inferior a la del reloj maestro, c, evitando las dificultades 
mencionadas. 

Tras N A1 ciclos cb t equivalentes a K-N A1 periodos T Xl el contador alcanza el 
valor terminal que alza la serial tca-1. En el siguiente ciclo de reloj cb se activa 
el flip-flop FF tipo RS (83), colocando la senal aj=1 . Este valor se propaga por 

20 el registro de desplazamiento SHRG (84) cuya longitud la establece el valor 
programado para N A2t accionado por el reloj maestro c. Transcurridos N A2 + 1 
periodos de reloj maestro, el valor aj=1 alcanza la salida del registro de 
desplazamiento, con lo que la salida pasa a HFD=7, habilitando a los circuitos 
de focalizacion dinamica una vez transcurridos N A periodos de reloj maestro, 

25 como indica la Ecuacion (46). De este modo, solo el registro de desplazamiento 
SHRG (84) debe operara la elevada frecuencia del reloj maestro. 

Tal estado permanece estable hasta que se haya realizado un numero N F de 
correcciones focales, programado en el contador CNT-F (86). Cada correction 
focal es indicada por la serial de entrada cef, que proporcion la logica de control 
30 CTRL (90) de la Figura 7 cuando se toman mutiples muestras entre focos, o es 
el reloj de muestreo CK k si solo se adquiere serial en los focos (Figuras 5 y 6). 
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Cuando CNT-F alcanza la cuenta terminal, se alza la senal tcf=1, provocando 
la puesta a cero del flip-flop FF (84) por su entrada R Al realizarse la funcion 
AND (81), la senal HFD retorna a 0, finalizando el proceso de correccidn focal 
progresiva. 

5 No obstante, la adquisicion de senales puede continuar con periodo de 
muestreo constante e igual al nominal v/m hasta completar el numero de 
muestras a adquirir, N St programado en el contador CNT-S (88). El 
rebosamiento de este valor activa la senal FIN que termina por complete la 
adquisicion de senal en este canal. 

10 Es importante observar que, finalizado el proceso de correccion focal 
progresiva al hacerse HFD=0, la direccion F de la memoria de codigos de 
correccion focal (75) queda congelada en el ultimo valor alcanzado por el 
contador (76) de las figuras 5, 6 y 7, lo que posibilita iniciar automaticamente la 
siguiente adquisicion con los valores de correccion focal almacenados a 

15 continuation, haciendo innecesario programar un nuevo valor inicial de F. De 
este modo es posible operar con diversas estrategias, como por ejemplo: 

a) Generation de imagen sectorial, donde cada adquisici6n se realiza con un 
angulo de deflexion diferente. Los codigos de correccion focal que 
corresponden a los angulos de deflexion sucesivos se almacenan 
20 secuencialmente en memoria, evitando tener que programar la direccion inicial 
F de la memoria de cddigos de correccibn focal antes de cada adquisicibn. A 
diferencia de otras tecnicas, esta permite mantener el error de muestreo 
acotado y constante para todos los angulos de deflexion y a todas las 
distancias. 

25 b) Realization de barrido lineal u otras formas de adquisicion donde una 
subapertura activa compuesta por un subconjunto de los elementos del array, 
cambia en cada adquisicibn, eliminando ciertos elementos e introduciendo 
otros. Analogamente, los codigos de correccion focal que corresponden a las 
configuraciones sucesivas se almacenan secuencialmente en memoria, 

30 evitando tener que cargar la direccion inicial F de la memoria de codigos de 
correcci6n focal antes de cada adquisicibn. 
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c) Estrategias mixtas, en las que se adquiere senal de multiples direcciones 
utilizando los mismos codigos de correction focal y diferentes valores N A (k) 
para diferentes direcciones con los que se adquiere el primer foco. En este 
caso, el valor de la direction inicial debe restablecerse con anterioridad a la 
5 adquisicion. 

El dispositivo APD (60), Figura 4, realiza las funciones de apodizacion y 
apertura dinamica. La apodizacion, que consiste en multiplicar la senal recibida 
en cada canal por un numero real positivo A k comprendido entre 0 y 1, mejora 
las caracten'sticas del haz conformado y el contraste en la imagen al reducir los 
10 lobulos laterales. La apertura dinamica permite mantener una resolution lateral 
constante operando con un numero F# constante, que requiere que el tamano 
de la apertura activa aumente progresivamente en el tiempo. Ademas, la 
apertura dinamica facilita la aplicacidn de la tecnica de Correcci6n Focal 
Progresiva en las regiones mas proximas al transductor. 

15 La Figura 9 muestra un diagrama de bloques del dispositivo APD (60). Para la 
funcion de apodizacion utiliza un multiplicador MULT (61) que obtiene el 
producto de cada uno de los valores de la secuencia e k por el coeficiente A k de 
apodizacidn asignado al canal k. En una realizacibn preferente, este coeficiente 
es un numero entero sin signo de a bits. Si los valores e k se expresan con w 

20 bits, de los a+w bits que componen el producto se mantienen los w bits mas 
significativos con redondeo en funcion del valor del primer bit ignorado y el 
signo de e k , incluyendose la logica correspondiente en el propio dispositivo 
MULT (61) para obtener la salida h k . La funcion realizada es: 



25 Debe hacerse la observation de que la funcion de apodizacion puede 
habilitarse o no con sencillos circuitos de control no mostrados en la figura 9. 
Basicamente, la salida h k es la entrada e k si la funcion de apodizacidn esta 
inhabilitada, respondiendo a la Ecuacibn (47) en otro caso, funcibn que se 
realiza con un simple multiplexor gobernado por la senal de habilitacion de la 

30 apodizacion. 




JU 



(47) 
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Para la funcion de apertura dinamica, se dispone un contador binario CNZ (62) 
que se carga inicialmente con el valor MZ(k) o numero de muestras nulas 
expresado por la Ecuaci6n (25) para el canal k. Este contador se acciona con 
un reloj cuya frecuencia es la nominal de muestreo, cs. Su operacion se habilita 
5 con la activacion de la senal HFD=1 que inicia el proceso de focalizacion 
dinamica y es proporcionada por el dispositivo HFOC (80), Figuras 4 y 8. Por 
su parte, la senal HFD=1 actua sobre la senal S del flip-flop FF (64) poniendo la 
senal z=7 que selecciona la entrada nula (63) al multiplexor MUX (65). 

En este estado, se inicia la adquisicion y procesamiento de senales en todos 
10 los canales. En aquellos en los que el contador CNZ (62) no haya alcanzado la 
cuenta terminal (senal fc), las muestras de salida en f k seran de valor nulo, por 
lo que el sistema opera como si el elemento asociado no estuviera activo. Tras 
MZ ciclos del reloj de muestreo, durante los que se han producido otras tantas 
muestras nulas sobre f kt se alza la senal tc-1 que pone a cero el flip-flop FF 
15 (64) y habilita la salida de los valores h k obtenidos por el multiplicador (61). De 
este modo se proporciona un numero programable de muestras nulas al 
comienzo de cada adquisicitin en cada canal, lo que permite programar las 
caracteristicas de la apertura dinamica. La secuencia f k es: 

f k = {0, 0, ...,0, A^NZ+1), AifiKfNZ+2), ...} (45) 

20 donde e«(i) es la muestra / de la secuencia e k y que, por consiguiente, se 
corresponde con la muestra / de cualquier otra secuencia. De este modo se 
mantiene el orden para realizar posteriormente la composicion coherente de 
todas las sefiales, sean o no nulas. 

25 Otro de los aspectos importantes de la presente invencidn es su capacidad de 
realizacidn modular, cuyo concepto se muestra en la Figura 10. Una realizaci6n 
preferente agrupa en un modulo (50) todas las funciones de procesamiento de 
senal para un subconjunto de PS elementos del array, asi como otras 
funciones de control global. En particular, un modulo incluye la memoria MJ 

30 (59) que proporciona los codigos de alejamiento focal J, junto con los circuitos 
necesarios para su lectura y carga inicial. Ademas, cada modulo (50) contiene 
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S submodulos (40), donde cada submodulo process las senales de P 
elementos, y un combinador (55) que permite enlazar multiples modulos en 
cadena. 

Internamente, los S subm6dulos (40) que conforman un mbdulo, se encadenan 
5 para realizar la composition coherente de los PS canales en cascada a traves 
de los FIFOs (52). Esta estructura facilita la operation en paralelo de todos los 
canales de procesamiento digital, donde los FIFOs compensan 
automSticamente las diferencias en los instantes de adquisicion de senal en los 
diferentes canales. 

10 Para ello, cada uno de los elementos (11) del array (10) proporciona la senal 
recibida a uno de los P procesadores elementales, PE, (30) integrados en un 
submodulo (40). Cada procesador elemental PE realiza el acondicionamiento 
de senal y digitalizacion en el instante de llegada de la senal al elemento que le 
corresponde, siguiendo la metodologia descrita anteriormente, e incluye los 

15 dispositivos que fueron mostrados con detalle en las Figuras 4 a 8. 

Cada PE proporciona a su salida una secuencia ordenada de valores 
correspondientes a las senales recibidas por el elemento asociado desde cada 
uno de los focos o muestras intercaladas. Estas salidas se suman (41) para 
realizar la composici6n coherente de las P senales que corresponden a un 
20 submodulo. 

Por su parte, el dispositive (42) suma los resultados obtenidos localmente en el 
submodulo con los que le haya colocado el precedente en su cola FIFO (52). 
De este modo, a la salida de cada submodulo se obtiene una secuencia 
ordenada que corresponde a las senales focalizadas de sus P elementos mas 
25 los de los subm6dulos que le preceden en la cadena, en el orden de llegada de 
senales desde cada foco. El primer submodulo de la cadena recibe un valor '0' 
(51) en sustitucion del FIFO (52), considerando que no existe un submodulo 
precedente. 

A su vez, el combinador (55) suma en (58) los resultados obtenidos por este 
30 modulo y que se almacenan en el FIFO (56) con los que el modulo anterior ha 
colocado en la cola FIFO (57). Como en ambos casos los resultados son 
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secuencias ordenadas de valores correspondientes a focos sucesivos, la 
secuencia de salida del sumador refleja este mismo orden. El encadenamiento 
de L mtidulos (50), cada uno de ellos con S subm6dulos (40) con P 
procesadores elementales en cada submodulo proporciona un sistema de 
5 focalizacion dinamica con la tecnica de Correction Focal Distribuida para 
N=LSP elementos. 

Es, por tanto, una estructura que presenta diversas ventajas: 

• Es modular, facilitando el desarrollo de circuitos y mayor flexibilidad de 
aplicacion. 

10 • Procesa en paralelo, lo que permite alcanzar elevadas tasas de calculo. 

• Tiene una capacidad ilimitada de expansi6n desde el punto de vista l6gico. 

La Figura 11 muestra una de las posibles realizaciones del sumador (41) 
integrado en un submodulo (40). Este sumador recibe las P secuencias ej a e P 
correspondientes a los valores de las muestras adquiridas por los P 
15 conversores AID pertenecientes al sub-modulo. En el caso mostrado en la 
figura es P=8, siendo la estructura generalizable a cualquier otro valor, 
preferiblemente potencia de 2. 

Las muestras que corresponden a cada foco en las secuencias e 1 a e P se 
obtienen en forma ordenada, esto es, la primera muestra obtenida por los P 

20 conversores A/D corresponde al primer foco, la siguiente al segundo foco y asi 
sucesivamente. Este razonamiento es extensivo aplicado a las muestras 
obtenidas entre focos, con o sin la opci6n de Correcci6n Focal Progresiva 
Variable. Sin embargo, los instantes de muestreo no son simultaneos en todos 
los canales, sino que es posible que un canal adquiera varias muestras 

25 sucesivas mientras que otros esperan a que las senales correspondientes 
lleguen al elemento asociado. 

Este desfase se corrige autom£ticamente con las memorias FIFOs distribuidas 
(43), (45) y (47) junto con los sumadores (44), (46) y (48). El FIFO (49) realiza 
la misma funcion pero con referenda a las secuencias de muestras de dos 
30 submodulos. 
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Los FIFOs se asocian en parejas, cuyas salidas se suman cuando ambos 
disponen de, al menos, una muestra almacenada. Esta situacion se detecta 
con un sencillo circuito de control que ejecuta una lectura simultenea en ambos 
FIFOs cuando no estan vacios y escribe el resultado obtenido por el sumador 
5 correspondiente en el FIFO conectado a su salida. Cada lectura libera una 
posicion de memoria en los FIFOs, que puede ser ocupada con la escritura de 
una nueva muestra. Cada escritura ocupa una nueva posici6n de memoria en 
el FIFO correspondiente. Las lecturas en los FIFOs asociados son simultaneas, 
pero no asi las escrituras que dependen de los instantes de muestreo en los 
10 respectivos conversores A/D, a su vez derivados del instante de llegada de la 
senal desde cada foco al elemento. 

La capacidad H de cada FIFO (43) debe ser suficiente para almacenar las 
muestras adquiridas hasta que se disponga de una muestra en el canal 
asociado, en cuyo momento opera el sumador (44), liberando la posicion 
15 ocupada por las muestras correspondientes y transfiriendo el resultado de la 
suma al FIFO (45). 

Analogamente, los FIFOs (45) deben tener una capacidad 2-H, los FIFOs (47) 
una capacidad 4 H, y asi sucesivamente. De este modo, la cantidad de 
memoria distribuida en forma de FIFOs requerida por un submodulo con P 
20 entradas es: 

HP + 2-HPI 2 + 4'if * PI 4 + .. . = H P \og 2 P (49) 

En el caso en que se utilizara un unico FIFO por cada canal, su capacidad 
deberia ser de H (P-1) posiciones, con lo que el total de memoria requerida 
seria HP-(P-1). Asi, por ejemplo, si H=76y P-8 t segun el metodo preferido 
25 cuya estructura se muestra en la Figura 6 se requieren 384 posiciones de 
memoria distribuida, mientras que en el caso de utilizar un unico FIFO por cada 
canal sen'an necesarias 896 posiciones. 

El FIFO (49), por su parte, almacena transitoriamente los resultados obtenidos 
por el submodulo actual hasta que el sumador (42) de la Figura 10 pueda 
30 liberar una posicion cuando el FIFO (52) disponga de un dato correspondiente, 
entregando en tat caso el resultado al submodulo siguiente. 



WO 2005/073957 



PCT/ES2005/070008 



41 

Breve description del contenido de las figuras. 

La Figura 1 muestra la geometria utilizada como ejempio, con un array 
compuesto por N elementos numerados 1, 2, 3, ... k, N-1 y A/, sobre el eje x. 
Los focos se situan a intervalos AR a partir de una distancia minima Ro. Se 
5 muestra el segmento entre un foco situado en coordenadas polares (R, Q) y el 
elemento k del array situado en coordenadas cartesianas (x k , 0), cuya longitud 
es L k (R,Q). 

La Figura 2 muestra la estructura general de un sistema de conformacibn digital 
de haces en reception generico. Se compone de un array (10) con N 
10 elementos (11) y un canal de procesamiento de senal por cada elemento. Cada 
canal contiene circuitos de amplificaci6n y acondicionamiento de senal AAS 
(12) y un digitalizador A/D (13). Las salidas de los N digital izadores son 
procesadas por el conformador CONF (14). 

La Figura 3 muestra la estructura utilizada frecuentemente para realizar la 
15 composition coherente de senales segun metodos conocidos, Consta de un 
sistema de retardo grueso (22) y otro fino (23), cuya operaci6n se gestiona 
desde un controlador local (24) que utiliza el contenido de una memoria de 
focalizacion dinamica (25). Las N salidas retardadas fi, f2, se suman en 
(20) para obtener la secuencia resultante r. 

20 La Figura 4 muestra la estructura de un canal de procesamiento digital 
preferentemente utilizado por el sistema objeto de esta patente. Contiene el 
elemento de amplificaci6n y acondicionamiento de senal AAS (12), el 
digitalizador A/D (13), una memoria FIFO (43) que almacena transitoriamente 
las muestras adquiridas, un generador de reloj de muestreo (70) que constituye 

25 un elemento fundamental de la presente invencion y un modulo de habitation 
de la focalizacion dinamica y control, HFOC (80). 

La Figura 5 muestra la estructura de principio para realizar el generador del 
reloj de muestreo (70) con codigos de correction focal Q expresados con 1 bit 
en la memoria MEM (75). El contador CNT (76) establece la direction de 
30 acceso a esta memoria. El reloj de muestreo CK kl propio del canal /c, se genera 
mediante logica distribuida en el registro REG-A (71a), el registro de 
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desplazamiento de longitud programable SHR-A (72), el flip-flop FFA (73), el 
multiplexor MUX (74) y el flip-flop FFB (77). 

La Figura 6 muestra la estructura de principio para realizar el generador del 
reloj de muestreo (70) con cbdigos de correccibn focal Q expresados con 
5 multiples bits en la memoria MEM (75). El contador CNT (76) establece la 
direccion de acceso a esta memoria. El reloj de muestreo CK*, propio del canal 
k, se genera mediante logica distribuida en el registro REG-A (71a), el registro 
de desplazamiento de longitud programable SHRA (72), el registro de 
desplazamiento de longitud programable SHRB (78) y el conjunto de inversores 
10 (75a). 

La Figura 7 muestra la estructura de principio para realizar el generador del 
reloj de muestreo (70) con codigos de correction focal Q expresados con uno o 
mas bits en la memoria MEM (75). El contador CNT (76) establece la direccion 
de acceso a esta memoria. En este caso puede insertarse un numero fijo o 

15 variable de muestras entre focos, segiin el contenido del registro-contador 
RCM (71b) controlando la variable ng que establece la longitud del registro de 
desplazamiento SHRB (78). La variable ng, de 1 bit, es proporcionada por el 
control CTRL (90) mostrado con m£s detalle en la parte inferior de la figura. 
Esta unidad consta de un contador de modulo programable CMOD (91), cuyo 

20 contenido y el valor del codigo de correccion focal Q, forman la direccion a una 
memoria local MDM (92). La salida g de esta memoria es invertida por (93) 
para obtener la serial ng. En este caso, existe la option de proporcionar el 
codigo de alejamiento focal J desde la misma memoria MEM de codigos de 
correccion focal (75), como se muestra con trazos de puntos o, 

25 preferiblemente, desde el exterior y a traves de una memoria FIFO (79). Por lo 
dem£s, el reloj de muestreo CK k del canal k se produce mediante la Ibgica 
compuesta por el registro REG-A (71a), el registro de desplazamiento SHRA 
(72) y el registro de desplazamiento SHRB (78). La unidad de control CTRL 
proporciona, ademas, la salida cef utilizada para habilitar el conteo de focos en 

30 la unidad HFOC (80) mostrada en las Figuras 4 y 8. 

La Figura 8 muestra esquematicamente el principio en que se basa la unidad 
HFOC (80) que realiza el control general de la adquisicion. Consta de un 
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contadorCNT-A (82) que se carga mediante Ida con el valor N A i> se habilita 
con la senal cea y se acciona con el reloj de baja frecuencia cb. La cuenta 
terminal tea actua sobre la entrada S del flip-flop FF (83). Su salida a/se 
retrasa una cantidad de ciclos de reloj maestro c expresada por N A 2 mediante 
5 el registro de desplazamiento SHRG (84). La puerta AND (81) produce un nivel 
alto en la salida HFD cuando sus dos entradas estan en alto. El flip-flop (83) se 
repone a 0 cuando se alza su entrada R al activarse la serial tct 

La serial tcf es la cuenta terminal del contador CNT-F (86) que puede 
inicializarse con un numero de focos Nf determinado al activarse la senal Idf, se 
10 habilita con la serial cefy se acciona con el reloj de baja frecuencia cb. Por su 
parte, el contador CNT-S (88) cuenta las muestras adquiridas, activando su 
salida FIN cuando se alcanza el numero programado A/ s . Esta programacion se 
realiza activando la senal Ids, el contador se habilita con ces y se acciona con 
el reloj de baja frecuencia cb. 

15 La Figura 9 muestra esquematicamente el principio en que se basa la unidad 
APD (60) que realiza las funciones de apodizacion y apertura dinamica. Consta 
de un multiplicador (61) que obtiene el producto de la entrada por el del 
coeficiente A k y entrega su resultado h k a un multiplexor MUX (65). Asimismo, 
dispone de un contador CNZ accionado por un reloj cs y habilitado por la senal 

20 HFD cuyo contenido inicial es /tfz, esto es, el numero de muestras nulas a 
aportar por este canal para realizar la funcion de apertura dinamica. Al alzarse 
la serial HFD se pone en alto la serial z de salida del flip-flop FF (64), 
seleccionandose la salida '0' entregada por la unidad (63) por el multiplexor 
MUX (65), que la proporciona como salida f k . Cuando se alcanza la cuenta 

25 terminal tc se activa la entrada R del flip-flop FF (64), pasando a nivel bajo su 
salida z con lo que el multiplexor MUX (65) selecciona la salida ft* del 
multiplicador MULT (61), entregandola por la salida f k . 

La Figura 10 muestra esquematicamente la arquitectura modular del sistema 
objeto de esta patente. Un modulo (50) contiene multiples submodulos (40), 
30 cada uno de ellos con varios procesadores elementales PE (30) cuyas entradas 
son las senales recibidas por los elementos (11) del array (10). Las salidas de 
los procesadores elementales son sumadas en (41) y la salida de este sumador 
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se utiliza como uno de los sumandos de una cadena de sumadores (42) y 
FIFOs (52). El primer subm6dulo sustituye este FIFO por una entrada nula (51). 
Los resultados proporcionados por el sumador (42) del ultimo submddulo se 
almacenan transitoriamente en un FIFO (56) para ser, a su vez, sumados en 
5 (58) con los resultados de un modulo precedente que han sido almacenados en 
el FIFO (57) en el sumador (58). Adicionalmente, el modulo contiene una 
memoria MJ (59) que proporciona los c6digos de alejamiento focal J a los 
diferentes procesadores elementales PE (30). 

La Figura 11 muestra esquematicamente el principio en que se basa una 
10 realization preferida del sumador (41), con un ejemplo para 8 canales, donde 
multiples FIFOs (43) almacenan transitoriamente las muestras obtenidas en 
cada procesador elemental e k . Las salidas de los FIFOs (43) son, a su vez, 
sumadas por pares en los sumadores (44), almacenandose el resultado en los 
FIFOs (45). Las salidas de estos FIFOs son sumadas en (46) y el resultado se 
15 almacena en los FIFOs (47). Analogamente, las salidas de estos FIFOs son 
sumadas en (48), almacenandose los resultados finales en el FIFO (49), cuya 
salida r representa la suma de las muestras e* a e 8 entregadas por los 8 
procesadores elementales que constituyen un submodulo en este ejemplo. 

La Figura 12 muestra un esquema de un ejemplo realizado de un modulo (200) 
20 del sistema objeto de esta patente. Dispone de conectores (120a) y (120b) por 
los que se introducen las senates analogicas que son distribuidas por los buses 
(121a) y (121b) a los conversores A/D dobles (101). Las muestras 
proporcionadas por estos conversores son introducidas por (106) en una FPGA 

(102) que realiza las funciones de procesamiento digital correspondientes a un 
25 submodulo (100), entre las que destacan la generation de los relojes de 

muestreo (105) y el suministro de resultados al siguiente submddulo por el bus 

(103) . El ultimo submodulo en la cadena suministra los resultados obtenidos a 
una FPGA (111) por el bus (104) para ser combinados con los resultados 
procedentes de otro modulo por el bus (113), suministrandose los datos 

30 resultantes por los buses (112) o (114). Estos buses esten presentes en el 
conector (110). A su vez, la FPGA (111) entrega senales de control (115a) y 
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(115b) a los circuitos de amplification y acondicionamiento analogico a traves 
de los conectores (130a) y (130b), respectivamente. 

La Figura 13 muestra las medidas experimentales realizadas sobre el modulo 
descrito en la Figura 12, mostrando el error absolute en ns para cada uno de 
5 los 32 canales. 

EJEMPLO DE REALIZACI6N DE LA INVENCI6N: 

En todos los casos descritos se obtiene el reloj de muestreo de cada canal 
utilizando muy pocos recursos hardware, fundamentalmente registros, 
10 contadores y una memoria para almacenar los codigos de correction focal, 
realizandose la adquisicion de las senates con una alta resolution temporal. 

Debido a la alta eficiencia de la metodologia descrita en cuanto a utilization de 
recursos hardware y, a diferencia de otras realizaciones, las tecnicas CFP y 
CFPV permiten integrar multiples canales de procesamiento en un unico 
1 5 circuito VLSI, junto con otras funciones descritas, incluso en dispositivos de uso 
general como son las FPGAs. 

Asi, en una realization demostrativa de esta invention, un modulo contiene 4 
subm6dulos, cada uno de ellos con 8 procesadores elementales. Este pequeno 
sistema puede operar con una apertura activa de hasta 32 elementos en 
20 barrido lineal o sectorial con apertura dispersa y arrays de 128 elementos. Con 
4 modulos se configura un sistema que opera con aperturas densas o 
aperturas activas de hasta 128 elementos. El sistema es escalable, anadiendo 
32 elementos activos por cada nuevo modulo instalado. 

La Figura (12) muestra gr£ficamente la configuraci6n del mbdulo, realizado en 
25 una tarjeta de circuito impreso de 220x100 mm. En la figura se han omitido los 
circuitos correspondientes al preprocesamiento analogico de senates, con 
posible excitation de los transductores, asi como los de posprocesamiento que, 
en esta realization particular, se basan en un computador integrado al que se 
conecta una Unidad de Control e Interfaz (UCI) con el sistema descrito. 
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El modulo (200) recibe las senales analogicas de 32 canales proporcionadas 
por los circuitos de amplificacion y acondicionamiento por los conectores (120a) 
y (120b). Contiene cuatro submodulos (100) y un combinador realizado en la 
FPGA (111). El combinador recibe senales procesadas por un modulo anterior 
5 en la cadena por el bus (113) y envia los resultados al que le sucede por el bus 
(114). Los parametros se programan desde la UCI a traves del bus de control 
(112) por el que, tambi6n, el ultimo mddulo de la cadena entrega los resultados 
finales. Asimismo, la FPGA (111) entrega un conjunto de senales digitales por 
los conectores (130a) y (130b) que son utilizadas por la unidad de 
10 procesamiento analogico. Por lo demas, esta estructura refleja fielmente la 
arquitectura presentada en la Figura 10. 

Las senales anal6gicas introducidas por los conectores (120a) y (120b) se 
distribuyen a los submodulos mediante los buses (121a) y (121b), 
respectivamente, donde son digitalizadas por los conversores A/D (101). La 
15 salidas digitales de los conversores A/D se introducen en las FPGAs (102) 
mediante los buses (106). Los relojes de muestreo de dichos conversores (105) 
se generan en las FPGAs (102) con la tecnica de Correcci6n Focal Progresiva 
descrita en esta memoria. 

Por su parte, los submodulos se conectan entre si mediante los buses 
20 bidireccionales (103). Junto al bus (104) que conecta el ultimo submodulo al 
combinador (111), estos buses se utilizan en un sentido para la programacion 
de parametros desde el de control (112) y en el sentido opuesto para 
suministrar resultados a los submodulos posteriores en la cadena o al 
combinador (111). 

25 El combinador realiza la suma coherente de los resultados obtenidos en el 
modulo actual a traves del bus (104) con los que recibe de un modulo anterior 
por el bus (113), suministrando los resultados de este proceso por el bus (114) 
al modulo siguiente o por el bus (112) a la Unidad de Control e Interfaz. Una u 
otra opcion son programables. 

30 En esta realization se han utilizado 16 conversores dobles de 10 bits, tipo 
AD9218 de Analog Devices para los dispositivos (101), 4 FPGAs tipo XC2S200 
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de Xilinx Inc. para los dispositivos (102) y una XC2S150 de Xilinx Inc. para el 
dispositivo (111). 

El reloj del sistema es de 40 MHz, velocidad a la que operan smcronamente las 
comunicaciones entre submodulos, mtidulos y con la UCI. Internamente, parte 
5 de los circuitos tambien operan a esta frecuencia o una fraccion que determina 
la de muestreo. Sin embargo, en las FPGAs (102) se multiplica por 4 la 
frecuencia del reloj del sistema para obtener el reloj maestro de 160 MHz. La 
resolucibn temporal es, pues, T x = 6.25 ns. y el error de muestreo est£ 
te6ricamente acotado a I e I < 3.725 ns. Con esto, pueden procesarse senales 
10 analogicas de hasta 10 MHz con \x > 16 y hasta 20 MHz con \i > 8, siendo [i la 
relation entre el periodo de la senal y la resolution temporal que afecta al 
rango dinamico y contraste de las imagenes. 

Cada FPGA (102) contiene los circuitos necesarios para realizar todo el 
procesamiento digital de 8 canales, incluyendo: 

1 5 • Apodizacion, con coeficientes 0<A k < 255/256 o deshabilitada (A k = 1, V k). 

• Apertura dinamica, con 0 < MZ k < 2 12 . 

• Deflexion del haz con foco fijo, con 0 < NA k < 2 14 . 

• Focalizacion dinamica segiin la tecnica de Correccion Focal Progresiva 
expuesta, con 4096 codigos de correccion focal de 1 bit por cada canal. 

20 • Programaci6n del numero de muestras a adquirir, con 2<N S < 2 17 . 

• Programacion del numero de muestras entre focos, con 4 < m < 64. 

• Programacibn del numero de focos por adquisici6n, con 0<N F < 2 12 . 

• Programacion de la frecuencia de muestreo de salida, con 2.5 < fs < 40 
MHz. 

25 • Circuitos de gestion de ahorro energetico y habilitacion de funciones. 

• Circuitos de programacibn de par£metros y procesamiento de resultados. 

• Circuitos para la verificacion del modulo, integrando un generador digital de 
senal. 
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Por su parte, el combinador contenido en la FPGA (111) tiene circuitos para 
realizar: 

• La gestion del modulo, transfiriendo valores de parametros a los 
5 submodulos por (104) y (103). 

• La recopilacion de resultados procedentes del modulo anterior por (113). 

• La composici6n de estos resultados con los del propio mbdulo, recibidos por 
(104). 

• El suministro de resultados combinados al m6dulo siguiente por (114). 
10 • El suministro de resultados finales a la UCI a traves del bus (112). 

• Circuitos para compensar el offset de las senales. 

• Circuitos para controlar la unidad de pre-procesamiento anal6gico. 

La Figura (13) muestra medidas experimentales relativas al error absoluto en 
15 ns. (ordenadas) del instante de muestreo en cada uno de los 32 canales (en 
abcisas) para un foco situado a una distancia de 50 mm., con un angulo de 
deflexion de 15°, velocidad de propagation c= 5900 m s 1 y array lineal de 5 
MHz con elementos uniformemente distribuidos a intervalos d=0.8 mm. 

El maximo error teorico es de 3.125 ns., que se verifica en 28 de los 32 
20 canales. El mayor error en los 4 canales que superan el valor teorico se debe a 
diferencias en el trazado de pistas en el circuito impreso y en las FPGAs de 
este modelo experimental y, en cualquier caso, no es significativo. El error 
cuadratico medio, mas determinate de la calidad de la imagen, es de 7.9 ns. 

Es evidente que pueden introducirse mejoras a la realization descrita puesto 
25 que existen dispositivos con mayor capacidad de integraci6n. Asimismo, 
pueden realizarse otras modificaciones que afectan al tipo de conversion A/D, 
distribution de dispositivos, etc., ya que el modelo descrito no es mas que un 
ejemplo demostrativo, sin que tales modificaciones supongan un cambio 
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sustancial a los fundamentos de la metodologia y sistema que han sido 
descritos en la presente memoria. 

Para la aplicacion del metodo, considerese un caso de imagen ultrasonica 
sectorial en pulso-eco, con un transductor lineal de N elementos centrado en el 
5 origen, con un tamano de aperture D. Se pretende cubrir con la imagen un 
sector de -45°<Q<45°. 

Si es la frecuencia de la sefial recibida, para verificar con suficiente margen 
el criterio de Nyquist, se elige una frecuencia de muestreo nominal f s = 4-f R . 
Ademas, para mantener un elevado rango dinamico y contraste en la imagen, 
1 0 se elige = fy/f R - 16, siendo f x la frecuencia del reloj maestro. 

Considerando el trayecto en ida y vuelta del pulso ultras6nico, el intervalo entre 
muestras es: 

AR=cT s /2 (50) 
Sustituyendo el valor AR dado por la Ecuacibn 1 resulta, con m = 1\ 

15 V=^ = 4 (51) 

Cuando se insertan m muestras entre focos, el valor de v = 4 m. Aplicando 
ahora la Ecuaci6n 18, con \x k \ < D/2 t donde D es el tamano de la apertura, la 
distancia minima para aplicar la T6cnica de Correccibn Focal Progresiva con 
apertura completa es: 



20 



D 2 



mv -2a 

cosG 



2-yJamv - a 1 



+ ^|sen9| (52) 



donde a=2^-7, siendo b el numero de bits con los que se expresan los codigos 
de correccion focal, con la condicion mv > 2a. 

La Tabla II muestra los resultados de la Ecuacion (52) para b-1, 2 y 3 en 
funcion del numero de muestras (fijo) entre focos. 

25 Debe observarse que el valor R(/D representa el mi'nimo valor del numero F# 
con el que se puede operar utilizando la funcion de apertura dinamica. En 
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general, se acepta como valido un valor F# m/n < 1 , que se obtiene con la 
mayoria de las combinaciones. 

Tabla II 

Relacibn FUnvn = R(/D para aplicar la TCFP con apertura completa {-45°< 9 < 
5 45°) 



m 


b = 1 


b = 2 


b = 3 


1 


0.6 






2 


0.8 


0.4 




4 


1 


0.6 


0.4 


8 


1.4 


0.9 


0.6 
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REIVINDICACIONES 

Esta patente descubre el metodo y sistema para conformation de haces en 
recepcion con una nueva tecnica denominada de Correccion Focal Progresiva, 
con la opci6n Variable. Segun dichos m6todo y sistema, se realiza la 
5 composition coherente de senates procedentes de fuentes situadas a lo largo 
de una direccion seleccionada, que Ilegan a un conjunto de A/transductores 
receptores o elementos distribuidos espacialmente, fonmando un array. Un 
conformador de haces en recepcion realiza una composicion coherente de las 
senales recibidas que proceden de fuentes puntuales situadas a lo largo de 

10 dicha direccion. La composicion coherente de las senales originadas en cada 
punto y recibidas por cada uno de los N elementos del array proporciona una 
unica serial de salida que representa a la que se obtendna con un transductor 
receptor unico, del mismo tamano que el array, focalizado individualmente en 
todos y cada uno de los focos, tecnica que se conoce en el ambito como 

1 5 focalizacion dinamica con deflexion. 

Las senales recibidas son amplificadas y acondicionadas analogicamente antes 
de ser digitalizadas. La tecnica de correccion focal progresiva, con o sin la 
opci6n variable, y el sistema que la realiza descritos en esta patente, se aplican 
a las senales digitalizadas por un conversor analogico-digital. Se denomina 
20 canal al conjunto de dispositivos que procesan la senal recibida por un 
elemento del array hasta su combination con las correspondientes a otros 
elementos. 

A diferencia de otros m6todos que forman parte del estado del arte, el metodo y 
sistema objeto de esta patente no se basa en retrasar las senales recibidas por 

25 cada uno de los elementos del array, ni utiliza procesos de interpolation. Antes 
bien, se fundamenta en muestrear las senales recibidas en los instantes en que 
Ilegan a cada uno de los elementos del array desde cada uno de los focos. El 
m6todo garantiza que la composicibn de muestras es coherente, esto es, s6lo 
se combinan las que corresponden al mismo foco y, ademas, no nedundante, 

30 adquiriendo unicamente las que van a utilizarse en el proceso de composicion. 
Requiere un reloj de muestreo de las senales no uniforme y diferente para cada 
uno de los canales. 
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Por otra parte y, a diferencia de otras tecnicas en las que cada reloj de 
muestreo queda establecido por una secuencia de T y '0' en una memoria 
individual para cada canal o se calcula en tiempo real mediante circuitos 
digitales especificos a partir de los valores de una diversidad de parametros 
5 programados durante la adquisicion, con el metodo y sistema objeto de esta 
patente los instantes de muestreo se determinan interpretando el contenido de 
tablas asociadas a cada canal, que se calculan y almacenan en memorias 
antes del inicio de la adquisicion de serial, haciendo innecesaria la carga de 
multiples parametros en tiempo de adquisicion. Ademas, la cantidad de 

10 almacenamiento requerida es limitada al indicar, en forma codificada, la 
correction focal que debe aplicarse al reloj de muestreo para obtener el 
siguiente foco. M£s aun, cada cbdigo de correccidn focal puede ser compartido 
por multiples muestras consecutivas. Finalmente, la arquitectura de la tecnica 
aqui patentada presenta una elevada regularidad, facilitando su integration en 

15 dispositivos VLSI y su realization modular. 

Se reclaman, como propias, las siguientes reivindicaciones: 

1. Un metodo de conformation de haces en reception, denominado Tecnica 
de Correcci6n Focal Progresiva, para aplicaciones en las que una 
diversidad de N transductores, agrupados en un array, reciben energfa 

20 vibratoria desde focos F it /=0, 1, 2,..., situados a lo largo de una direction 
programable, caracterizado porque muestrea la senal recibida en cada 
canal en el instante de llegada desde cada foco con un error absolute 
inferior a medio periodo de un reloj maestro de periodo T Xi de forma que las 
muestras consecutivas e*/, /=0, 1, 2, ... obtenidas en cada canal k 

25 correspondan a los valores consecutivos de la senal originada en cada foco 
F b i=0, 1, 2,..., con lo que cada muestra de la senal resultante r h /=0, 1, 2, ... 
se realiza mediante la suma de las N muestras e^, k-1 , 2, .., N : 

N 

^ i = 0,l 5 2,... 

o de una combinacion lineal como: 
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15 



N 



k=\ 

donde los coeficientes A k son valores fijos o programables, operation que 
se denomina suma o composition coherente de las N senates, 

2. Un metodo de conformation de haces en reception segun revindication 1 
5 caracterizado, ademas, porque para el canal receptor correspondiente al 
elemento k del array, que incluye circuitos para el acondicionamiento y la 
amplification de serial, el instante de muestreo que corresponde al primer 
foco F 0 se determina, en tiempo no real, a partir del numero de ciclos N A (k) 
del reloj maestro que median entre un origen de tiempos elegido con la 
10 restricci6n de que todos los tiempos sean positivos y el instante de llegada 
de la senal del foco Fo al canal /c, de modo que si el tiempo de vuelo de la 
serial correspondiente al foco F\ al elemento kes Ti(k), calculado a partir de 
la geometria del sistema y de la velocidad de propagation de las ondas en 
el medio, el valor NA(k) viene dado por: 



NAk) = 



T 



donde [-]u representa la funcion de redondeo al entero mas proximo, que 
determina N A (k) con un error absoluto inferior a Tx/2, funcion que puede 
realizarse mediante un contador de N A (k) ciclos de reloj maestro, 
generandose un pulso de adquisicion de senal en el canal k en el momento 
20 en que dicho contador rebosa; tal y como se ha descrito en la memoria, 
existen diversas formas de realization fisica, prefiriendose una en la que el 
contador se desdobla en dos partes, una de pequeno tamano accionada 
por el reloj maestro y la principal por un reloj cuya frecuencia es un 
submultiplo de la del reloj maestro 

25 3. Un m6todo de conformaci6n de haces en reception segun reivindicaciones 
anteriores caracterizado, ademas, porque para otros focos F b i > 0, el 
instante de muestreo en el canal k se determina en tiempo real a partir del 
valor Q*/de b de bits, almacenado en una memoria de correction focal 
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asociada a dicho canal k que codifica el adelanto a aplicar al valor nominal v 
de periodos del reloj maestro que hay entre dos focos consecutivos, donde 
el codigo Q# se calcula en tiempo no real a partir de la diferencia AT^ entre 
el instante de llegada al elemento kde las senales procedentes del foco i-1 
y el foco /como: 



e«=v- 



T 



donde [-]u representa la funcion de redondeo al entero mas proximo y v es 
el intervalo nominal entre focos, expresado en periodos del reloj maestro, 
programandose consecutivamente los codigos Qw que corresponden a 

10 focos sucesivos F h M, 2, 3, ... en una memoria de correction focal 
asociada al canal k cuyo contenido se lee en tiempo de adquisicion de 
serial para producir el adelanto variable en el instante de muestreo para 
cada foco F/ en funcion del valor de Qw mediante circuitos digitales de los 
que en la memoria descriptiva se han dado varios ejemplos, basados en 

15 registros de desplazamiento de longitud programable, siendo posibles otras 
alternativas, como el uso de contadores programables, sin que por ello se 
modifique la idea base expuesta en esta revindication. 

4. Un metodo y sistema de conformation de haces en reception segun 
reivindicaciones anteriores caracterizado porque permite insertar un numero 

20 fijo o programable m de muestras entre focos, por medio de un contador 
programable m6dulo N M =m que se incrementa cada vez que se adquiere 
una nueva muestra y, cuando retorna a cero, se incrementa la direccibn de 
acceso a la memoria de correction focal, facilitando que un total de m 
muestras compartan el mismo codigo de correction focal Qw, lo que permite 

25 aumentar la eficiencia de la codification que puede ser de una fraction de 
bit por muestra adquirida, manteniendo el error del instante de muestreo 
basicamente acotado dentro de medio periodo del reloj maestro al 
disponerse de circuitos digitales que distribuyen uniformemente el adelanto 
Qid entre las m muestras que se adquieren entre focos, circuitos que se 

30 basan en una tabla en la que se programan, antes de la adquisici6n y para 
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cada posible valor Q w < 2*-*/, un total de m bits g que indican, para cada una 
de las m muestras, si debe adelantarse (g=7) o no {g=0) en 1 ciclo de reloj 
maestro el instante de adquisicion de dicha muestra, constituyendo un caso 
particular de la metodologi'a expuesta en las reivindicaciones anteriores. 

5 5. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque puede variarse la distancia 
entre focos mediante la adquisicion de un numero variable de muestras m 
entre focos que se incrementa a partir de un valor inicial mo cuando el valor 
de un codigo de alejamiento focal, J h de 1 bit, asociado al foco F/es J, = 1, 

10 opcion denominada tecnica de Correccion Focal Progresiva Variable, que 
aumenta la eficiencia de uso de la memoria de correccidn focal al aumentar 
con el tiempo la distancia entre focos, siendo el valor de Ji comun a todos 
los canales para cada valor de /, pudiendo tabularse de forma totalmente 
centralizada para todos los canales, parcialmente centralizada afectando a 

15 un conjunto de canales, como se realiza preferentemente, o individualmente 
en cada canal, lo que permite establecer un balance entre la flexibilidad de 
aplicacion y la eficiencia de memoria; para su realization basta sustituir el 
registro del valor constante m asociado a cada canal k por un contador que 
se carga con el valor inicial m 0 antes de iniciar cada adquisicidn y, en 

20 tiempo de adquisicion, dicho contador se incrementa cada vez que se 
reciben senates desde un nuevo foco F, y el codigo J/=7. 

6. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque dispone de contadores 
para determinar el numero de focos N F y de muestras N s a adquirir, 

25 teniendo en cuenta que cada foco puede incluir la obtencion de m muestras, 
siendo m constante o variable, de forma que cuando el contador N F rebosa, 
se bloquean los mecanismos de la correccion focal progresiva y de 
actualizacion de los codigos de alejamiento focal, congelandose, ademas, la 
direcci6n de memoria que accede a los c6digos de correcci6n focal en cada 

30 canal, aunque la adquisicion de senates puede continuar hasta que el 
contador de muestras a adquirir N s alcance el valor terminal, facilitando que 
una unica memoria por cada canal pueda almacenar los codigos de 
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correction focal para multiples angulos de deflexion, multiples 
configuraciones de la apertura activa o una combinacion de ambas, 
pudiendo programarse en cualquier caso la direccibn inicial de la memoria 
de correccion focal antes de cada adquisicion u obtenerse de forma 
5 automatica si la secuencia de adquisiciones se corresponde con el orden en 
que se han almacenado los codigos de correccion focal para la actual 
aplicacibn. 

7. Un metodo y sistema de conformaci6n de haces en recepcion segiin 
reivindicaciones anteriores caracterizado porque opcionalmente puede 

10 operar con apertura dinamica, activando un contador asociado a cada canal 
/c, en el que se programa un numero N z (k) de muestras de salida de valor 
nulo que se producen a partir del comienzo de la adquisicion para su 
combinacibn con las muestras correspondientes de los restantes canales, 
de modo que, para estas N z (k) muestras, el resultado equivale a eliminar el 

15 elemento k en la composicion de las correspondientes muestras de salida, 
realizando efectivamente la funcion de apertura dinamica. 

8. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado, adem£s, porque puede insertar 
una funcion de apodizacion o multiplication de los valores de cada muestra 

20 adquirida por una constante programable para cada canal. 

9. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque realiza la composicion 
coherente mediante una combinacion de colas o memorias FIFO 
distribuidas y sumadores, preferentemente formando una estructura en 

25 arbol en la que cada dos FIFOs se asocian a un sumador de 2 entradas. 

10. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque facilita su realization 
modular con una arquitectura flexible que permite la composicion coherente 
de un numero indefinido de senales desde el punto de vista l6gico, para lo 

30 cual, cierto numero de canales se agrupan para formar un sub-modulo, un 
numero de sub-modulos configuran un modulo y un sistema se configura 
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con cierto numera de modules, donde todos los canales son iguales y 
realizan el procesamiento en paralelo de la information adquirida. 

11. Un metodo y sistema de conformation de haces en reception segun 
reivindicaciones anteriores, que utiliza la metodologia y modelos 
5 constructivos expuestos en la memoria y figuras adjuntas. 
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